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Abstract X-radiation from energetic electrons is the prime diagnostic of flare-acceler- 
ated electrons. The observed X-ray flux (and polarization state) is fundamentally a 
convolution of the cross-section for the hard X-ray emission process(es) in question 
with the electron distribution function, which is in turn a function of energy, direc- 
tion, spatial location and time. To address the problems of particle propagation and 
acceleration one needs to infer as much information as possible on this electron distri- 
bution function, through a deconvolution of this fundamental relationship. This review 
presents recent progress toward this goal using spectroscopic, imaging and polarization 
measurements, primarily from the Reuven Ramaty High Energy Solar Spectroscopic 
Imager {RHESSI). Previous conclusions regarding the energy, angular (pitch angle) 
and spatial distributions of energetic electrons in solar flares are critically reviewed. 
We discuss the role and the observational evidence of several radiation processes: free- 
free electron-ion, free-free electron-electron, free-bound electron-ion bremsstrahlung, 
photoelectric absorption and Compton back-scatter (albedo) , using both spectroscopic 
and imaging techniques. This unprecedented quality of data allows for the first time 
inference of the angular distributions of the X-ray-emitting electrons using albedo, im- 
proved model-independent inference of electron energy spectra and emission measures 
of thermal plasma. Moreover, imaging spectroscopy has revealed hitherto unknown 
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details of solar flare morphology and detailed spectroscopy of coronal, footpoint and 
extended sources in flaring regions. Additional attempts to measure hard X-ray polar- 
ization were not sufficient to put constraints on the degree of anisotropy of electrons, 
but point to the importance of obtaining good quality polarization data. 

Keywords Sun: flares; Sun: X-rays; Sun: acceleration; Sun: energetic particles 
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1 INTRODUCTION 

X-ray emission, because it is produced promptly in an optically-thin environment, 
is one of the most direct methods with which to study energetic electrons in solar 
flares. Such remote radiation measurements are generally functions of photon energy 
e, direction fl, and time t. A common description of incoherent and partially polar- 
ized X-ray radiation typical of solar flares is in terms of its photon intensity I{e,fl,t) 
(photons cm~^ sr' keV"^ CGS units are used throughout the paper) and the 
fractional degree and orientation of linear polarization (P, tf'). The emergent photon 
intensity is the number of photons dN that escape from the source in the time interval 
t to t + dt in the energy range e to e -I- de, from a solar source with a direction and 
angular size defined by the elementary cone dfl, oriented with respect to the normal 
of the detector of area dS: 

dN = I{€,n,t)d€dSdndt. (1.1) 

In other words, the intensity is the number of photons emitted by a unit solid angle 
source (sr) on the Sun and received in a unit of time (1 s) in a unit of energy band 
(1 keV) by a unit detector area (1 cm^) at the Earth. 

For optically-thin radiation, the emergent intensity of the radiation is simply a lin- 
ear convolution of the cross-section for the pertinent emission process and the electron 
phase-space distribution function /e(v,r, t) (electrons cm^"^ [cm s^^]^"^), or equiva- 
lently the particle flux F{E,fl',r,t) (electrons cm^'^ s~^ sr ~^ keV~^), differential 
in energy E, and velocity solid angle fl' . Since dE = mevdv in the non-relativistic 
regime, it follows that F{E,fl',r,t) = fe{y,r,t)/me, where me (g) is the electron 
mass. 

For an elementary bremsstrahlung source of ambient plasma density n(r), located 
at position r on the Sun, along the line of sight 17, which is subjected to an electron 
flux spectrum F{E,fl',r,t), the emergent photon flux spectrum at distance R is the 
convolution 

I{e,fl,t)=[ I I n(r)F(£,r,n',t)Q(n,fl',e,£;)d£'dn'd^, (1.2) 
J I Jn' Jt 

where I is the distance along the line of sight, Q{Q,,d' ,e, E) (cm^ keV^^ sr^^) is 
the cross-section for the pertinent hard X-ray emission process(es), differential in e 
and n. The dominant hard X-ray emission process in solar flares is bremsstrahlung 
radiation associated with electron deceleration in the Coulomb field of an ion or other 
electron. For this process, the angular dependence of Q depends only on the angle 
6' — n'fJ between the incoming electron f2' and the emitted photon f2 directions, so 
that Q = Q{e,E,e'). 

To deduce the physical properties of energetic particles in solar flares from observed 
hard X-ray quantities, the electron flux spectrum F{E,r,fl',t) must be deconvolved 
from the emission cross-section Q{e, E, 6') in this integral. Although the relation be- 
tween the observable quantity I{e,fl,t) and the physical electron flux F{E,r,fl',t) is 
linear, it is still nontrivial and the deconvolution requires some rather insightful tech- 
niques. This fundamental problem and the progress toward the solution of the problem 
using spectroscopic, imaging and polarization measu rements from th e Reuven Ramaty 
High Energy Solar Spectroscopic Imager (RHESSI. iLin et al.ll20o3 ) are addressed in 
this review. 
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In Section [J] we review the physical processes leading to X-ray emission in so- 
lar flares, including free-free electron-ion bremsstrahlung (Section 12. ip . free-bound 
electron- ion emission (Section 12. 2p . and free-free electron-electron bremsstrahlung 
(Section 12. 3p . This review, however, does not discuss bound-bound transitions, nor 
the emissions due to energetic ions such as nuclear gamma-ray lines, ion gamma-ray 
continuum, and pseudo-continuum. 

Since the emission mechanisms are well established, hard X-rays (HXR) are often 
viewed as one of the most direct (i.e., least affected by propagation effects) diagnostics 
of solar flare electrons. However, it must not be forgotten that downward-propagating 
X-ray photons are effectively scattered toward the observer by electrons in the dense 
layers of the solar atmosphere, thereby complicating the diagnostic potential of hard 
X-ray radiation. Section [3] reviews spectroscopic and imaging techniques to infer this 
Compton backscattered (albedo) component. Progress in the deduction of primary, 
i.e., directly flare-emitted, and photospherically-reflected X-rays from the observed 
spectrum, is discussed. 

Section 13] presents an overview of the results obtained using purely spectroscopic 
data (i.e., data integrated over the source volume). Both forward fitting and regularized 
inversion techniques to deduce the energy dependence of the mean source electron flux 
spectrum F{E) (averaged over volume and solid angle) are reviewed. Properties of 
the electron flux distribution, deviations from power-law forms, low-energy cutoffs and 
interpretation in terms of a thermal source are discussed. 

Section [5] provides methods and results for the case of an anisotropic angular dis- 
tribution F{E,ft) of electrons. The role of X-ray Compton scattering in the solar 
atmosphere (solar albedo) in deducing this angular distribution is discussed. Recent 
RHESSI polarization measurements and their implications for electron anisotropy are 
also reviewed. 

Section [6] focuses on the spatial structure F{E,r) of the electron flux, using imag- 
ing spectroscopy observations from RHESSI. Recently-developed, visibility-based tech- 
niques to optimize the inference of electron maps are discussed. Section [7] highlights 
the major flnding from RHESSI and discusses the open questions. 

The role of these re sults in the co n text of multi-wavelength observations of so- 
lar flares is discussed bv lFletcher et al The implicatio ns of these find i ngs fo r 
electron transport, and a cceleration models ar e discussed by iHolman et al.l (|201ll '). 
iHannah et all (|201ll '). and lZharkova et all l|201ll '). 



2 X-Ray emission processes and energetic electrons 

When energetic electrons are deflected in close encounters with ambient particles (both 
electrons and ions), a bremsstrahlung (literally, "braking radiation") photon is pro- 
duced. As stated in Section [T] the cross-section Q{e,E,d') for electron-ion free-free 
bremsstrahlung is a function of the emitted photon energy e, the pre-coUision elec- 
tron energy E, and the angle 6' between the directi on of the pre-coUision electron and 
the outgoing photon (see, e.g.. iKoch fc Mot3ll959l ). For simplicity we often consider 
only the solid-angle-integrated form of the cross-section Q(e, E) (and the corresponding 
scalar electron flux F[E, r, t]); however, it must be remarked that the effects of the an- 
gular dependence of the cross-section can, for highly-beamed electron distributions, be 
quite significant, leading to substantial differences in the number of electrons required 
to produce a given hard X-ray fiux - see Section [5] 
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The photons of energy from a few keV to a few hundred keV under considera- 
tion are mostly produced by coUisional electron-ion bremsstrahlung in the solar at- 
mosphere. Bremsstrahlung emission from energetic electrons is more efficient than 
inverse Compton scatt ering or synchrotron emission from the same electron pop- 
ulation ( Korcha5ll967l ). The responsible electrons have kinetic energies _E up to a 



few hundreds of keV, and so will coUisionally stop within a column density A'^ ~ 
E'^ ~ 10^^ < 10 ^^ cm~^, wh ere e (esu) is the electronic charge and A the 
Coulomb logarithm (e.g., lEmsliel Il978l '). Thus, the column density in a solar coro- 
nal loop of density 10^'^ cm~^ and length 10^ cm is 10^^ cm^^, which stops elec- 
trons up to ~10 keV, while the upper chromosphere with typical column densities 
10" -10 cm can stop electrons with the energies of 10 — 300 keV. The cross- 
section for scattering the emitted photons is of the order of the Compton cross-section 
ac « nro « 2 X 10~ cm , where ro = e /meC is the classical electron radius, so that 
the optical depth t = N ~ 10""* and there is negligible self-absorption in the source, 
i.e., the source is optically thin. In Section [31 we shall address the issue of so-called 
"albedo" photons - photons that are emitted downward toward the solar photosphere, 
which is optically thick, and subsequently backscattered toward the observer. For now, 
we consider only "primary" photons, i.e., those initially directed toward the observer. 
Integrating the photon intensity J(e, 12) over the solid angle subtended by the source 
of area A, dQ = dA/R^ and making the volume element substitution dfl dl = d^r/R^, 
the fundamental Equation H1.2|) shows that the observed bremsstrahlung flux (at time 
i; hereafter understood) at the Earth (photons s~^ keV~^ cm~^ of detector area) is 

/(e) = J^n{r)F{E,r)Q{e,E)dEd\ (2.1) 

where the second integral is taken over the source volume V . 



2.1 Electron-ion bremsstrahlung 

For electron- ion br emsstrahlung, the fu ll form of the cross-section Q{e,E) is given 
by formula 3BN of iKoch fc Motj (|l959l ). N umerical com putations may be facilitated 
by use of the simplified form published bv iHau j ( 199?! ). A frequently-used analytic 



approximation to the cross-section Q{e,E) is the Kramers form 

Q{e,E)^Z^^, (2.2) 

where ao = (8q/3) (meC^) = 7.9 x 10"^^ cm^ keV. Here a ~ 1/137 is the fine 
structure constant, me is the electron mass, Z is the mean ion charge, c is the speed of 
light. A more accurate analytic form, valid in the non-relativistic limit, is the Bethe- 
Heitler form 



Q(e,E) = Z^^ln : ' " (2.3) 



ao, 1 + yi - e/E 



eE 



l-^Jl^T/E 



For purely spectral observations, a spatially-integra ted form of the basic Equa- 
tion (|2.H) is appropriate. In this case, we can write (e.g.. lBrownlll97ll l 

[nVF{E)\Q{e,E)dE, (2.4) 
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where n = (1/V) fyn{r)d^r and F(E) = (1/nV) fy n{r) F{E,r) (fir. The quan- 
tity F(E) (electrons cm~^ s^"'^ keV~^) is termed the mean electron flux spectrum 



(IBrown et al. 2003); it has also been termed the "X-ray emitting electron spectrum" 



bv l Johns fc LinI (|l992aH bh. Since the quantity n V (the number of target particles in the 



emitting volume) is dimensionless, the units of the quantity \nV F{E)] are the same as 
those for electron flux, viz., electrons cm^^ s^^ keV""'^. For large events, typical values 
of F(E) and \nV F{E)] at a representative energy _E ~ 20 keV are of order 10^® and 
10^'^, respectively. 

The inference of [n V F{E)] corresponding to an observed /(e) may be accomplished 
in several ways, which are described in detail later in this chapter. A review of the 
different degrees of effectiveness of these techniques in discerning the overall magnitude 
of, overall spect r al sha pe of, and form of "local" features in F{E) has been presented 
bv lBrown etall l|2006l '). 



It is of crucial importance to note that the quantity [n V F{E)] is the only quantity 
that can be inferred unambiguously (i.e., without additional model assumptions) from 
the source- integrated bremsstrahlung emission /(e). Use of Equation (|2.4p to obtain 
\nV F{E)] for a given /(e) is therefore a fundamental issue in the interpretation of 
solar hard X-ray spectra. Oncejn V F{E)] has been determined, the actual magnitude 
of F{E) depends on the valued of n and V. 

As an example of the use of the mean source electron spectrum to determine phys- 
ical properties of solar flares, let us consider the inference of the accelerated electron 
flux spectrum J-'q{Eq). The bremsstrahlung yield, the number of bremsstrahlung pho- 
tons emitted between e and e + de from an electron of initial energy Eq in a plasma of 
density n(r), may be written 

/■^'' nir)Q{e,E)viE)dE 

where dE/dt is the energy loss rate (here assumed a function of E only). For energy 
losses in a cold target due to binary collisions with the background electrons (radiation 
energy losses are much smaller and ca n be ignored) , dE/dt = ~{K/E) n(r) v{E), where 
K = 2ne'^ A = 2.6x 10"^* cm^ keV^ (|Emslielll97i ). and so the total observed flux from 



an injected distribution with energy spectrum J^q^Eq) (electrons cm ^ s ^ keV ^) is 



POO 



/(e) = Wo) H^, Eo) dEo 

fEQ 



47r/i2 K 



- / MEo)dEo / EQ{e,E)dE, 



(2.6) 



where A (cm ) is the area of the flare. Reversing the order of integration in (|2.6|l gives 

^ . POO POO 

^W = T-R2l7/ EQ{e,E)dE To{Eo)dEo, (2.7) 

and comparing this with the fundamental Equation 1)2. 4|) yields the result 



^ The astute reader will note that as the source volume V — > oo, the value of F(E) — > 0. 
This formal difficulty may be removed in practice by the truncation of the emission volume V 
at some reasonable upper limit. 
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— A 1"^ 

nVF{E) = -E To{Eo)dEo. (2.8) 

Jeo=E 

From this it follows straightforwardly that 

permitting the determination of the quantity ATq{Eq) (electrons s~"'^ keV~^); this 
qua ntity represents the rate of injection of electrons per unit energy (e.g..l Holman et"al] 
12003! ). Clearly J'q{Eq) is a nonnegative function, and so Equation (|2.9[) constrains F{E) 
to either be a decreasing function of E or, at worst, a function that increases more slowly 
than E. As we shall see below, certain recovered forms of \nV F{E)] fe.g.. IPiana et all 
can have difficulty satisfying this constraint; however it is also shown that such 
difficulties m ay be removed if th e effects of photospherically-backscattered (albedo) 
photons (e.g.. iKontar et al]|2006l ) are taken into account in Equation (|2.4|) . 



2.2 Free-bound emission 

In fitting or inferring mean source electron flux spectra F{E), free-bound recombi- 
nation emission by nonthermal electrons had alw ays been neglec ted c ompared with 
free- fr ee electron- ion bremsstrahlung, as argued by iKorchakI ( 1967h and lLandini et all 
(|l973l ). For hot plasma hard X-ray sources (coronal or in soft X-ray footpoints) this 
is inconsistent with inclusion of recombination as si gnificant for thermal electrons of 
similar energies in thermal spectrum modeling (e.g.. ICulhane|[l969l : ICulhane fc ActonI 
\m(i ). Of importance also is the fact that the estim ated coronal abundance Az for Fe 
is now much higher than in (e.g. iLandini et aPl 19731 ) . The recombination emission rate 
oc Az for hydrog enic ions of charge Ze and abundance Az with AzZ'^ « 1 for H 
and ~ 40 for Fe^^"*". [Srown et al.l ( 2010l ) have therefore re-examined the importance of 



nonthermal electron recombination. 

In the hydrogenic Kramers approximation, the free-bound emission rate from a 
plasma of proton density rip and volume V from nontherm al electrons with mean 
source electron spectrum F(E) is given by (jBrown et al.ll20ld ) 



where x ~ 13.6 eV is the H ionization potential, Z^^g the effective charge on the ion, 
re the classical electron radius and a the fine-structure constant, n is the principal 
quantum number of the empty shell into which the electron recombines and pn is a 
'vacancy factor' which takes into account the ratio of available to total states in that 
shell. [Note that Equation (pTTOll applies for e> Ec + Z^^ xir? ; for e < -I- ^cff xl'n^ , 
Jr{() = 0.] Since the recombination cross-section falls off as , it is adequate to take 
n — n„iin, the value for the first empty n-shell. A key feature of free-bound emission, 
unlike free-free, is that for any specific shell, each electron energy value E maps to a 
unique photon energy value e so that J_R(e) is a much more direct reflection of F{E) 
than is the bremsstrahlung Jb{^) convolution of F{E), the former preserving features 
in F{E). 
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Electron distribution 




Ptioton cJistribution 



Electron distribution 




Fig. 2.1 Left: the thin-target (including both thermal and nonthermal) electron and photon 
(right) spectra for two different plasma temperatures, 20 MK (top) and 30 MK (bottom). 
Electron nonthermal spectra for various values of low energy cutoff Ec and electron spectral 
index 5: Ec = 1 keV and <5 = 2 (blue and dashed); Ec = 10 keV and 5 = 5 (black and solid); 
The corresponding photon spectra: blue solid line for 5 = 2, Ec = 1 keV with free-bound 
emission; green dashed line for 5 = 2, Ec = 1 keV without free-bound emission; black solid 
line for 5 = 5, -Ec = 10 keV with free-bound; red dashed line for 5 = 5, -Ec = 10 keV without 
free-bound. 



Comparing Equa tion (I2.10p with the corresponding Kramers free-free expression, 
iBrown &: MaUikI (|200& . ,2009) concluded that nonthermal free-bound emission is neg- 
ligible in cold sources. However, they concluded that in hot plasmas (T ~ 10 — 30 
MK), such as coronal sources and the soft X-ray component of footpoint sources, free- 
bound emission can dominate for steep spectra and low cut-of f energies (Figure [2.ip . 
Subsequently, a serious error in iBrown fc MaUikI l|2008l . [ 2OO9I ). involving the use of 
incorrect values of ^ in their Equation (13) (i.e., Equation |2.10) above) was rec- 
ognized ( Brown et al. 20ld ). This substantially reduced both the magnitude of the 



cross-section and the energy shift e — E involved in recombination. Because of the 
typically steep electron spectra involved, the latter substantially affects the number of 
electrons responsible for emitting the photon in question. In the amended results Jfl.(e) 
is never dominant, even for for hot sources, but can account for up to ~ 30% of the 
flux J_B(e) + Jfl,(e) in the range ~10-30 keV for the non-thermal component dominating 
thermal as shown in Figure ETTI For such hot sources, the free-bound emission could be 
comparable with such effects as albedo and differences in bremsstrahlung cross-section 
used, and could be used to diagnose the sharp features in F(E) from J(e). Even more 
important is the fact that J_R(e) adds edges to the total J(e) with the result that, in 
data with good signal to noise ratio, inversion (essentially differentiation) of data J(e) 
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to yield F{E) based on bremsstrahlung alone (essentially differentiation of J) could 
result in spurious features in F{E) just as happens when albedo is ignored. 



2.3 Electron-electron bremsstrahlung 

Energetic electrons propagating in the solar atmosphere encounter ions and electrons 
(both free and bound) and hence can produce X-ray emission via both electron-ion and 
electron-electron bremsstrahlung. 

When the maximum electron energy is much larger than the photon energies under 
consideration, the photon spectrum resulting from a power-law sp ectrum of e lectrons 
F{E) oc E^^ is also close to the power-law form I{e) oc e~'^ { Hand Il989l ). How- 



ever, while for pure electron-ion bremsstrahlung 7 ~ 5 -I- 1, for pure electron-electron 
bremsstrahlung a significantly shallower photon spectrum, with 7 ~ 5, results. Thus, 
the importance of the electron-electron bremsstrahlung contribution increases with 
photon energy and the enhanced emission per electron leads to a flattening of the pho- 
ton spectrum /(e) above ~300 keV produced by a given F{E) (|Hauelll975h or, equiv- 



alently, a steepening of the F{E) form required to produce a given /(e). Kontar et al.l 
(|2007l ) provide a discussion of the essential differences between electron-electron and 
electron-ion bremsstrahlung processes. We here provide a succinct summary of that 
discussion and note that the properties of electrons with the energies above ~400 keV 
are also crucial for ion diagnostics and radio emission. 



As is well known l|Koch fc Mot j 1 19591 ). the cross-section for electron-ion brems- 



strahlung scales as . Further, when considering electron-electron bremsstrahlung, 
the possible binding of target electrons to their host ions in a neutral or partially- 
ionized medium is not significant. Hence, in a quasi-neutral target of particles with 
atomic number Z, the bremsstrahlung cross-section per atom for emission of a photon 
of energy e by an electron of energy E is in general equal to 

Q{e, E) = Z2(3e-p(e, E) + ZQe-e{e, E), (2.11) 



where Qe-p{^,E) and Qe-e{^,E) are the cross-sections, in the lab oratory frame, 
for ele ctron -proto n , and electron-electron bremsstrahlung, as given by iKoch fc Mot j 
(II959I) and iHaud (|l989l ). respectivelj0. It is also important to note that while the 
electron-ion cross-section is finite for all e < E, the cross-section for electron-electron 
bremsstrahlung vanishes above a maximum phot on energy, due to the necessar ily fi- 
nite energy carried by the recoiling target electron dHaudl 19751 : iKontar et al ] l2007l l. For 
highly non-relativistic electrons {E meC^), emax — ^ E/2; only for hig hly relativistic 
electrons {E ^ me(?) can emax — >■ E (see Figure [2^ . As discussed bv lKontar et al.l 
( 20071 ). :his result has important implications for the form of the photon spectrum 



produced by electron-electron bremsstrahlung. 

In Section 21 we discuss the application to the recovery of the electron spectrum at 
pertinent (mildly-relativistic) energies. 



Note that ther e is a typographical error in the form of Qe-e in IHaud jligsgD : see 
IKontar et"all \200t ) for details. 



10 



E.P. Kontar et al. 




Fig. 2.2 Maximum photon energy emax produced by electron-electron bremsstrahlung, ex- 
pressed as a fraction of tlie incident electron energy E (in units of the electron rest mass mc^), 
for various values of 6, the angle between the incoming electron and the outgoing photon tra- 
jectories. For clarity, only curves for 9 = 0°, 60°, 120°, and 180° are labeled; the curves for 
e = 30°, 90° and 150° lie between these Cafter iKontar et al.ll2007l '). 



3 Primary and Compton backscattered X-rays 

As discussed in Section [2l the atmosphere above the X-ray bremsstrahlung-producing 
region is optically thin. However, the very dense lower photospheric layers below the pri- 
mary source are not. Consequently, photons emitted downwards are efficiently Compton 
backscattered by atomic electrons in the photosphere. The observed hard X-ray (HXR) 
flux spectrum from solar flares is therefore a combination of primary bremsstrahlung 
photons with a spectrally modified component from photospheric Compton backscat- 
ter of the downward-directed primary emission. This backscattered component can be 
significant, creating new features and/or distorting or masking features in the primary 
spectrum, and so substantially modifying key diagnostics such as the electron energy 
budget. 

Photons of energy above ~100 keV penetrate so deeply that they are lost t o the ob- 
serve r, while below ~10 keV they are mostly photoelectrically absorbed (e.g.. iTomblinI 
\l9Ti ]. Therefore, the reflectivity of the photosphere has a broad hump in the range 10- 
100 keV, with a maximum around 30-40 keV. At some energies and view angles the re- 
flectivity approaches 90%, so the observed spectrum may be very substantially affected 
by backscatter. This effect is well-known in s olar physics (and more generally in X-ray 
astronomy (see iMagdziarz fc Zdziarskilll995l ). and there have been several discussions 
of its influence on ob served X-ray spectra (e.g.. lTomblinlE972l : ISantangelp e t al. I 19731: 



Bai fc Ramaty| l978[). and on the electron spectra inferred from them (e.g.,ljohnsfcTinl 



1992bl : [Alexander fc Brownlbood : iKasparova et al.ll2005l : | Kontar et al 



2006 



Understanding and modeling backscatter (albedo) has beco me even more im por- 
tant with the advent of high quality X-ray spectra from RHESSI dLin et aLlbOoi ) with 
spectral resolution as high as ~1 keV, in combination with uncertainties as low as a few 
percent (for strong flares). Generally, contamination of the observed X-ray spectrum 
by reflected photons leads to a flattening of the spectrum and hence to an underesti- 
mation of the electron s pectral index if the contribution of backscattered photons is 
not taken into account ( Bai fc Ramatvl[l978h . Extrapolation to low electron energies 
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using such an underestimated spectral index leads to a substantial underestimation 
of the total electron energy in a flare. Indeed, a low-energy cutoff in an uncorrected 
electron spectrum is n ot even required if a "true" primary (albedo-corrected) electron 
spectrum can be used ( Kontar et al.|[200^ 'l. 



3.1 Spectroscopy of the photospheric albedo 

Downward-emitted photons are either absorbed or Compton-scattered, with some of 
the latter returned toward the observer, adding to the total flux detected (X-ray 
albedo). Scattering takes place on electrons, whether free or atomic. To account for 
elements heavier than hydrogen, the Co mpton cross-section is multiplied by an effec- 
tive mean atomic number Z = 1.2 (e.g.. lAschwandenll2005 l for clement abundances). 
The detailed density structure of the medium is irrelevant (Tomblin 1972). 

Absorption, on the other hand, does depend strongly on chemical composition, and 
the best estimate of photospheric abundances should be included. The heavy elements 
Fe/Ni play the most important role from 6 — 8 keV up to ~ 30 keV, while lighter 
elements contribute below 6 keV ( Morrison fc McCammonlll983l ). 



3.1.1 Green's junction approach 

Propagation, absorption, and Compton scattering of primar y hard X-rays can b e 
straightforwardly studied using Monte-Carlo simulations (e.g., iBai fc RamatvlFlOTSl ). 
This technique is ideal for obtaining the reflected, and hence the total observed, pho- 
ton spectrum for a given form of the primary X-ray spectrum. However, the primary 
spectrum is generally unknown and unlikely to be an exact power-law, as is sometimes 
assumed. Therefore , an approach independent of the primary spectrum is required 
( Kontar et aPbood l. 



For any isotropic primary spectrum /p(eo), (photons cm ^ s ^ keV ^), we can 
write the secondary, backscattered spectrum /5(e) as 



oc 

/ 



-^s(<:,m)= / -''p(<;o)G(^i, e, eo)rff;o, (3.1) 



where G(/i,e, eg) is an angle-dependent Green's function and /i = cos 6*, where 9 is 
the heliocentric position angle of the primary hard X-ray source, or the angle between 
the Sun center - observer and Sun center - X-ray source lines. The observed spec- 
trum, at photon energy e and direction p, is I pit) + Ip(e). The im portance of not 
averaging over viewing angle may be seen from lBai fc Ramatvl (Il978[). Using Green's 
functi o n analytical fi t s to M onte Carlo simulations derived bv lMaedziarz fc Zdziarskil 



(|l995l ). lKontar et al.l (120061 ) calculated functions G{p, e, eg) shown in Figure 13711 for so- 
lar flare parameters. Green's functions, as calculated, account for Compton scattering 
and bound-free absorption. 

The shape of the Green's function depends on the energy of the primary photon. For 
primary photons with low energies eq < 30 keV, the Green's function has a rather sim- 
ple structure close to a Dirac delta- function (Figure 13. ip . showing that backscattering 
is dominated by the first scattering (especially at low energies) , with the contributions 
from higher orders of scattering being generally small. 
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Fig. 3.1 Green's functions G{fi,e,eo) [ keV^^], including Compton scattering and photo- 
electric absorption, for three primary photon ener gies ep = 2 0, 80, 300 keV and fi = 0.7 
(6 45°) calculated usin g approximations dMaedziarz &: Zdziarski.,1995i ) for solar conditions 
(after iKontar et ani2006l'l . 



3.1.2 Reflected X-ray photon spectrum 

Because the reflectivity is spectrum-dependen t (Figure 13.211 . the albedo spectrum als o 
depends on the shape of the primary spectrum dKontar et al.ll200^ : lKotoku et al.ll2007l V 
While previous studies considered only the results for prescribed power-law or ther- 
mal p rimary spectra (cf . iBai fc Ramatvlfl978l : I Johns fc Linlll992bl : [Alexander fc Brownl 
I2OO2I ). the Green's function method allows more general forms of the primary spectrum. 
The total observed spectrum /(e) is given by 



I 



(3.2) 



For a measured I[e,^), we may obtain the primary spectrum Ip{e) by solving the 
integral equation ()3.2|l . In practice measurements yield discrete quantities and the 
integral equation 1)3. 2[) is used in the matrix form 



(3.3) 



where we have used the summation convention for repeated indices, and introduced 
the Green's matrix 



Gijifi) = J G{fi,ei,eo) deQ. (3.4) 

Due to sharp features in the Green's function (Figure [3.1[) . the integration in Equa- 
tion (|3.4p is best performed via a change of variable to the wavelength (reciprocal 
energy) domain. 
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Fig. 3.2 Upper panel: Primary (solid line), reflected (dotted line) and total (dashed line) 
photon spectra calculated assuming a primary spectrum Ip{t) oc , and using the Green's 
function for an X-ray source at heliocentric angle d = 45°. Lower panel: Reflectivity, defined 
as the ratio of reflected to primary fluxes R{e,d = 45°) = Iji{e) / Ip{t). The reflectivity taken 
from|Bai & Ramatv (1978) is shown with diamonds. Two a bsorption edges of Fc at 7.1 keV 
and Ni at 8.3 kcV are seen in the reflected component (after iKontar et al.ll2006l ). 



3.2 Imaging of photospheric albedo 

3. 2. 1 Expected spatial stgnatures of albedo 

Before considering observational approaches to the spatial isolation of albedo, it is 
instructive to consider an elementary model of the solar backscatteri ng process. To 
do this, we make four simplifying assumptions (e.g.. lBrown et allllOTSl ): first, that the 
primary bremsstrahlung is generated as a point source at a height h above a planar 
scattering surface; second, that this surface is perpendicular to the line of sight (e.g., 
disk center flare); third, that the primary X-ray emission and backscattering processes 
are isotropic; and fourth, that absorption of the scattered photons can be neglected. 
With such a scenario, the albedo source would extend out to the horizon as seen from 
the primary source. The albedo surface brightness would be determined by the distance 
from the primary source and the scattering location and by a cosine illumination factor 
and so the surface brightness would fall off as (1 -f [r//i]^)~'^/^, where r is the radial 
distance between the scattering location and the sub-source point on the scattering 
surface. Note that the scale of the radial profile is determined by the height of the 
primary source. 

The implications of easing some of the restrictions in this simple model are illus- 
trated in Figure 13.31 For a primary source with finite size or structure, the resulting 
albedo patch would be a convolution of the primary source with the aforementioned 
profile. For a primary source located away from disk center, the center of the albedo 
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Fig. 3.3 Left panel: Model of a single primary source with albedo contours. The solid contours 
at logarithmically-spaced intervals show the primary source down to 6.25% of the maximum, 
and the brightest portion of the albedo patch at 3.1% of the peak of the primary. The dashed red 
contours show the integrated flux of the albedo patch. Note that 50% of the albedo flux arises 
from an area about 10 times larger than the primary source, given its height of 12 Mm. Right 
panel: A model of a double source 10 Mm above the photosphere with the resultant albedo 
patch shown by flux contours (red d ashes). Ther e is considerable overlap of the primary source 
and albedo patch in this case (after ISchmahl fc Hurford 2002) . 



patch is displaced toward disk center and presents an elliptical shape oriented parallel 
to the limb. 

More riKorous calcul ation of the spatial properties of an albedo patch are provided 
bv lBai fc Ramatvl (|l978l ) who show, for example, that limb darkening would make the 



albedo more easily detected near disk cen ter. In addition, the al bedo patch will be 
energy- and primary-spectrum-dependent ( Kontar fc Jeffre"vll2O10l ). Nevertheless, the 



simplified model does suggest that the albedo has three potentially observable spatial 
signatures: first, for sufficiently high primary source altitudes, the albedo would be 
significantly larger in extent than the primary source, with a size scale that increases 
with source height; second the albedo source would be displaced toward disk center 
by a distance /isin^Q, where ^0 is the heliocentric angle; third, the albedo source 
would be elongated parallel to the limb with a minor to major axis ratio of cos 9q. 
Observationally, however, because of the relative size scales of the primary source and 
its albedo patch, the albedo surface brightness would be only a small fraction (typically 
only a few percent) of that of the primary source. This would pose a potential challenge 
for conventional imaging systems because of scattered light; for reconstructed RHESSI 
hard X-ray images, typically limited in dynamic range to about 10:1, this would seem 
to make the spatial detection of albedo even more problematic. 

3. 2. 2 The spatial-frequency signature of albedo 

The potential observational difficulties posed by the low surface brightness of the albedo 
source can be eased if we Fourier transform the source of size ~d with the distribution 
exp(— x^/d^) and consider the amplitude of the Fourier components as a function of 
spatial frequency oc exp(— fc^d^), where k is the spatial frequency. Compact primary 
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sources have Fourier amplitudes that fall off at high spatial frequencies kd ^ 1. For 
spatial periods large compared to the source dimension kd <^ 1, however, the ampli- 
tudes are effectively constant. Further, the Fourier transform of the (1 + [r//i]^)~^/^ 
profile of an albedo source implies an amplitude that varies as e ~^^ , where k is the 
spatia l frequency. As illustrated in the lower half of Figure 2 in ISchmahl fc HurfordI 
(|2002l ). this could be readily distinguished from that of the primary source. In effect, 
the Fourier transform integrates over the faint, distributed albedo source so that it 
becomes potentially detectable. The most obvious albedo signature would then be an 
"excess" in Fourier amplitude at low spatial frequencies over that expected for the 
compact primary source. The exce ss would be comparable to the reflected fraction, 



viz., up to several tens % (|Bai fc Ram atv 1978; Kontar ct al. 2 006 ). This is relevant to 



RHESSI, because with nine logarithmically-distributed spatial frequencies, its rotating 
modulation collimators (RMCs) directly measure the Fourier components of the source 
distribution. 

3.2.3 Detection of a "halo" component 

Applying such considerations to flare observations, early analysis of RHESSI data 
(fSchmahl fc Hurford ,2002 ) directly fit the observed modulated light curves to derive 
an average modulation amplitude for each RMC. This technique, although simple, 
assumed circularly symmetric sources and was subject to potential statistical issues 
since it directly fit the counts which were only sparsely populating the short time 
bins. Nevertheless, the analysis showed (Figure 13. 4|) that the modulation amplitude 
continued to increase tow ard the coarsest sub coUimator even though the latter had 
a 183" spatial resolution ( Hurford et al.l[2002l ). This was taken as evidence of a hard 



X-ray "halo" component, consist ent with expectations from albedo. 

The analysis was refined by ISchmahl fc HurfordI \2QQ± . who formed individual 



back-projectionimages with the nine RMCs. In this case, the peak calibrated intensity 
of each back-projection image corresponded to the azimuthally-averaged modulation 
amplitude for the corresponding spatial frequency. The results, illustrated in Figure [531 
confirmed the presence of a non-Gaussian large-scale component consistent with a 
halo source. This technique alleviates potential statistical concerns associated with 
directly fitting the sparsely-populated time bins, but still requires circular sources for 
unambiguous interpretation. 

3.2.4 Direct use of visibility measurements 

Since the initial reports of the halo component, analysis techniques have been developed 
that enable the RHESSI data to be directly transformed into calibrated measurements 
of the vis ibilities (specific Fourier components) as a funct ion of spatial period and ori- 
entation (|Schmahl fc HurfordI [iooi : ISchmahl et al.ll2007l ). There are several reasons 



why visibilities are more appropriate for albedo determination: the visibilities are fully 
calibrated so that instrumental issues can be cleanly separated from solar issues; the 
expected visibilities can be directly calculated from source models; each visibility mea- 
surement is independent and has well-determined statistical errors whose propagation 
can guide subsequent conclusions; visibilities can be readily calculated using code that 
is now an integral part of the RHESSI object-oriented software package; visibilities 
are determined from linear transforms of the observed count rates, so that they can 
be combined in time or energy as desired; visibilities are well-suited to more complex 
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Fig. 3.4 Relative amplitudes of 7 flares vs. spatial frequency in the energy range 12-25 keV. 
The departures of the profiles from a Gaussian shape are interpreted as a "halo" compo- 
nent. Since (for computational simplicity) the fits were assumed to be azimuthally symmet- 
ric, the profiles must be considered to be polar averages over 360 degrees of rotation (after 
ISchmahl fc Hurford|[200g) . 



sources since an observed visibility is the sum of the corresponding visibilities of its 
components. 

Figure 13.61 illustrates observations expressed in terms of visibilities. This event 
occurred at about 60° longitude, so the albedo emission is expected to be relatively 
weak, although there may be a signature of albedo for subcoUimators 8 and 9, as 
indicated by the excess observed flux over the model flux. 



3.2.5 Future prospects for vtsibiUty-based albedo measurements with RHESSI 

In the context of albedo, visibilities are used in two ways: first, to generate maps of the 
source using a direct visibility imaging method [such as back-projection (MEM-NJIT) 
( Schmahl et"ai]|2007l ) iiO obtain primary source positions and fluxes; and second, to use 
a visibility-based forward-fit algorithm to parameterize simple models of the primary 
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Fig. 3.5 A Fourier inversion produces the brightness profile in the energy range 12 — 25 keV. 
Since the albedo surface brightness is very low, the cumulative integrated flux from r to oo is 
plotted as a function of r. If the primary source profile is Gaussian, the integrated fiux profile is 
Gaussian. The clear departures from a Gaussian shape i n these profiles indicates t he presence 
of extended sources inferred to be albedo patches (after [Schmahl &: Hurfordll2003| ). 



source and the corresponding parameters of the albedo component (size, intensity, 
ellipticity, location), or maximum entropy 

The MEM-NJIT mapping algorithm is useful for obtaining the qualitative source 
configuration, although it is less reliable for quantitative measurement of source sizes. 
On the other hand the forward-fit algorithm requires a good starting point in parameter 
space, but can yield reliable source parameters. 

There are practical difficulties that remain to be overcome. As suggested by Fig- 
ure [3]6] the relative calibration (to within ~5 %) of the different subcoUimators (espe- 
cially the coarse ones) is critical to the isolation of albedo. Regrettably, the detector- 
to-detector calibration cannot currently support such an objective. However, using 
relative visibilities (viz., normalizing each RMC's visibilities to that detector's spa- 
tially integrated response) eliminates the effects of the detectors' relative efficiencies 
and so provides a potential approach to bypassing the calibration issue. Confirmation 
of this approach could be achieved with near-limb flares (presumed to be without sig- 
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Fig. 3.6 Observed visibility amplitudes (crosses vfith blue error bars) for a flare interval as 
a function of subcoUimator (SC=l-9) and position angle (PA = — 180°) of the grids in the 
energy range 12 — 25 keV. Each of the 9 vertical panels shows the amplitude as a function 
of PA for one subcoUimator (labeled by digits below the X-axis). The red curve represents a 
model using two Gaussian sources, and the squares show the residuals relative to the model. 
For a given subcoUimator (5 and 6 are good examples), the amplitude rises and falls while the 
grids rotate from PA=0 to PA= 180° . Such patterns can correspond to an extended or double 
source. 



nificant albedo) or at energies where the albedo component is minimal. Other potential 
improvements are corrections for azimuthal averaging in each visibility, improvements 
in the present forward fit algorithm, perhaps with the use of an alternative search 
algorithm and with the use of third harmonic to add additional spatial frequencies, 
visibilities 

The first flares for which albedo can be parameterized will necessarily be spatially 
simple (single, compact, and strong) and located within a few arc min of disk cen- 
ter. Eventually, it is expected that spectroscopic tools for albedo will be combined 
with imaging tools to provide comprehensive albedo information for a large subset of 
RHESSI flares. Since the albedo intensity and location depends on electron directiv- 
ity, the potential reward of spatially- based albedo diagnostics is well worth the effort 



a gn 

required to refine the analysis tools. I Kontar fc Brownl 2006bl ) 



4 The electron energy spectrum 

In this section, the angular dependence of the bremsstrahlung cross-section and the an- 
gular/spatial/temporal characteristics of the electron distribution are neglected. The 
primary source photon spectrum Ip{e) may be therefore treated simply as the convo- 
lution in electron energy of the solid-angle-averaged bremsstrahlung cross-section and 

the mean source electron flux (equation 12. 4|) . 

In general, we first correct for instrumental effects (see lSmith et al.l2002l : ISchwartz et al.l 

I2OO2I ) such as pulse pileufj^], then correct the observed hard X-ray spectrum for albedo 

^ Pulse pileup is an issue for large solar flare spectra with high count rate in RHESSI 
detectors. Un-physical counts are recorded when pairs (or more) of low-energy photons, arriving 
nearly simultaneously, are detected as a single energy count at higher energies (for details. 
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effects fSection 13. 1|) to obtain ttie primary source spectrum Ip{(-)- As expiained in 
Section [2l the functionai form of tlie primary plioton energy spectrum I{e) (subscript 
'P' hereafter understood) contains crucial information on the form of the mean source 
electron flux spectrum F{E), information that may in turn be used to reveal properties 
of the electron acceleration and propagation processes. 

We first discuss forward fitting fSection l4.1|) and regularized inversion fSection l4.2l) 
methods of extracting the electron energy spectrum F(E) from noisy hard X-ray data 
7(e). 



4.1 Forward fitting 

Forward fitting is the process of quantitatively comparing a parameterized model with 
observational data. Criteria are established to determine acceptable fits and, if the 
model is capable of providing acceptable fits, a best fit, that gives the most probable 
values of the model parameters, is determined by minimizing chi-squared (x^), the sum 
of the squares of the normalized residual^. 

Models may be based on the apparent structure of the data, a physical model, or a 
combination of the two. At the lowest X-ray energies, RHESSI flare spectra can con- 
tain thermal hremsstrahlung and free-hound (reco mbination) continua from the hottest 
plasma in flares, in addition to spectral lines fe.g.. IPhillips et aDbood '). The X-rays at 
higher energies (normally > 10-20 keV) are dominated by electron-ion bremsstrahlung 
from energetic, nonthermal electrons fSection l2.ip . Free-bound radiation from nonther- 
mal electrons may also contribute fSection 12. 2|l . At higher energies, electron-electron 
bremsstrahlung can become signiflcant fSection 12. 3|l . At 7-ray energies, spectral lines 
excited by energetic, nonthermal ions and positronium continuum emission can be 
present. The quality of a fit is also dependent on a careful subtraction of background 
counts before obtaining the spectral fit. Spectral fits are only reliable over the range of 
photon energies for which the flare emission is well above the background. 

RHESSI's high spectral resolution often makes it easy to distinguish the thermal 
component of an X-ray spectrum from the nonthermal component, especially in large 
flares. An example of this is shown in the top panel of Figure B7I1 a spatially-integrated 
spectrum from SOL2002-07-23T00:35 (X4.8); below a photon energy ~40 keV, the ther- 
mal component clearly dominates over the flatter nonthermal component. For many 
spectra, however, the thermal and nonthermal components are not so clearly distin- 
guishable. T his was the case du ring the early rise of SOL2002-07-23T00:35 (X4.8), 
for example ( Holman et ahlfioOSl ). For these spectra, the most likely model can some- 



times b e deduced from the time evolution of the flare spectra and/or flare images (cf. 
ISui et al. 2005). 

The thermal component of the spectrum can typically be well fltted with the form 
of /(e) from an isothermal plasma at temperature T : /(e) oc {EM/eT^^'^) exp(— e/fcT), 
where k is Boltzmann's constant and EM is the emission measure (J neUidV). The 



see lDatIowdll976bl lah. Current pileup corrections for spatially integrated spectra have limited 
precision that might be inadequate for the events with extremely high count rates. As of April 
2011, no standard image pileup corrections yet exist. 

* A normalized residual is defined as the difference between the measured and model- 
predicted value, divided by the uncertainty in the measured value. The division by the uncer- 
taint y gives greater wei ght to measured values with smaller relative uncertainties (e.g.. Chapter 
15 of lPress et aIlll992D 
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Fig. 4.1 Fit to a spatially-integrated spectrum from SOL2002-07-23T00:35 {X4.8). Top Panel: 
Photon flux spectrum (plus signs) integrated over the 20 s time period 00:30:00—00:30:20 UT. A 
fit to the spectrum {solid curve) consisting of the sum of the bremsstrahlung from an isothermal 
plasma {dotted curve) and the bremsstrahlung from a double power-law mean electron flux 
distribution with a low-energy cutoff. Middle Panel: Residuals from the flt in the top panel 
(observed flux minus model flux divided by the Icr uncertainty in the observed flux). Bottom 
Panel: Best fit mean electron flux distribution times the m ean plasma density and source 
volume, plot ted as a fu nction of electron energy in keV Cafter lHolman et al.|[200a) . Note that 
pulse pileup llSmith et al...2002i : .Kontar et al...2003) might be an issue for this flare. 

fit in Figure 14.11 for example, gave a temperature of 37 MK and an emission measure 
of 4.1 X 10"*^ cm~^. A fit to thermal bremsstrahlung alone is often adequate, but tliis 
does not account for the spectral line complexes at ~6.7 keV and ~8 keV or for re- 
combination radiatior[f]. Tlie thermal component of RHESSI spectra is now routinely 
fitted (included in to standard RHESSI software) using tiie latest version of Cliianti 
( Landi et allbood ). which incorporates all the emission mechanisms important at low 
energies. The multithermality of plasma and the corresponding emission measure dif- 
ferential in temperature are addressed in Section [4.61 

^ RHESSI spectroscopy in t he range of energies below ~ 10 keV is of ten complicated by non- 
diagonal instrument response llSmith et al.ll2002l : [Schwartz et al ] |2002l'l and some instrumental 
features llPhillips et al.ll2006f) 
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The nonthermal component of the spectra can usually be fitted adequately with 
either a single or a double power-law photon flux spectral model. Sometimes a third, 
flatter power-law component is included at low energies to simulate a low-energy cutoff 
in the electron distribution . Such fits are useful for examining the evolution of flare 
spectra with time. They do not, however, contain any direct physical information about 
the electrons responsible for the observed emission, ft is therefore more interesting to 
fit the photon spectra with the radiation from a model electron distribution, typically 
assumed to have the form of a double power law with a possible low- and/or high-energy 
cutoff. This form allows sharp breaks in the electron distribution (either a mean electron 
flux F{E); Equation |2.4| or, for a thick-target model, an injected electron distribution 
Tq{Eq); Equation [2.7] '). However, due to the filtering of the bremsstrahlung cross- 
section Q{e,E), such breaks are generally smoothed out in the corresponding photon 
spectrum I{e). 

The nonthermal part of the flare spectrum in Figure l47l1 is flt with the bremsstrahlung 
from a double power-law mean electron flux distribution with a low-energy cutoff: 

r 0; E <Ec 

F{E) ^ I AE-^^; Ec<E<Ei, (4.1) 

[ AeI"""^^ E-^^- Ef, < E. 

(bottom panel). The highest value of the low-energy cutoff Ec consistent with a good 
flt to the data was used; the value of Ec is not constrained below this value because 
of the dominance of thermal radiation. This fit therefore provides a lower limit to the 
energy in nonthermal electrons. The spectrum could not be acceptably fit with a single 
power law; note that the location of the break energy E^, is at a higher energy than 
the apparent location of the break in the photon spectrum; this is because all electrons 
with energies above a given photon energy contribute to the radiation at that photon 
energy. 

The photon spectrum residuals (using a sum of the isothermal and nonthermal 
models) are shown in the middle panel of Figure 14.11 Besides providing a reduced 
(X^ divided by the number of degrees of freedom in the fit) close to 1, the residuals 
from a good fit should be random and uncorrelated and have a near-normal distribution 
N{0, 1). For the event in question, the residuals do not exhibit this desired behaviour, 
with significant deviation at photon energies between 10 and 15 keV. This is not a 
well-understood issue and the explanati on could be either due to unaccounted non- 
diagonal response fsee lSmith et al.ll2002l) or the presenc e of a "superhot" component 



in the thermal continuum spectrum ( Caspi fc Linllioiol ). Practically, broad spectral 



"line" is often included in the model to account for this feature; with the inclusion of 
this ad hoc feature, it is generally possible above ~10 keV to obtain good fits to the 
RHESSI spectra without assuming the presence of any systematic uncertainty in the 
data above the level of Poisson noise. 

A comparison of a poor forward fit (left panels) to a good forward fit (right panels) 
is shown in Figure 14.21 Both fits are to the same spectrum and in both the model is 
bremsstrahlung from an isothermal plasma plus a double power-law photon spectrum. 
In the fit on the left, the break energy Ef, is fixed at a value ~21 keV, so the model 
is a single power law above this energy. In the flt on the right, the break energy is 
allowed to adjust to a value that gives the best fit to the data. The reduced for the 
flt on the right is 1.0; in the fit on the left it is 1.4, still consistent with an "acceptable" 
flt. However, the "long wavelength" oscillation in the residuals for the flt on the left 
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Photon Energy (keV) Photon Energy (l<eV) 

Fig. 4.2 Albedo-corrected RHESSI spectrum (crosses with error bars) at the hard X-ray peak 
(11:44:28-11:44;32 UT) of SOL2002-06-02T20:44 (Ml.O). The solid line shows the combined 
isothermal (dotted line) plus double power-law (dashed line) spectral fit. The spectral fit before 
albedo correction is overlaid (gray, solid line). Left panel: The nonthermal part of the spectrum 
is fitted to a double power-law model with a break energy at ~21 keV. Right panel: The same 
spectrum is fitted to a double power-law model with a break energy at ~40 keV. The reduced 
values of the fits in the left and right panels are 1.4 and 1.0 , respectively. The normalized 
residuals are plotted below each spectrum (after [Sui et al.ll2007fl . 



provides the most obvious clue that the model is not adequate; the overall shape of 
the spectrum is not consistent with the assumed model form. Adding more parameters 
(in this case a variable break energy E^,) results in a more acceptable distribution of 
residuals, as in the right panel (Figure [4. 2|l . 

The gray curves in Figure are fits to the data before albedo, assumed to be from 
isotropically-emitted photons, is taken into account (see Section tn]) . They demonstrate 
the significant impact that albedo can have on the inferred spectrum of the emitted 
radiation. 

An inadequately resolved issue with forward fitting is determining the uncertainty 
in the fit parameters and the resulting model function. Since, in general, the fit is not 
linear, the fit parameters are not independent and the uncertainties are not necessarily 
even symmetric around the best-fit values. In Figure IT2l for example, notice that the 
temperature of the isothermal component is adjusted to a high er value to comp ensate 
for the low break energy in the double power- law component I Sui et a"l]|2007^ . The 
uncertainty in the value of the low-energy cutoff to the mean electron flux fit function in 
Figure l47l1 is small in the positive (higher-energy) direction, but indefinitely large in the 
negative direction! As long as a good initial choice is made for the fit parameters, the 
process of obtaining the best fits is relatively quick. An efficient method for determining 
the uncertainties in the fit parameters and function is not in place, however. Bayesian 
Monte Carlo approaches to determining these uncertainties are robust but slow. A 
practical solution to this important issue is badly needed. 



4.2 Regularized inversion 

The unprecedented energy resolution of RHESSI hard X-ray spectra has introduced 
the need, perhaps for the first time in solar hard X-ray spectroscopy, to apply sophisti- 
cated mathematical tools for information retrieval in order to fully exploit the physical 
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significance of the recorded spectra. It is well established that the most effective math- 
ematical framework f or this problem lies in the theory of linear inverse problems (e.g., 
ICraig fc BrOTrolll986l '). In this setting, linear integral equations of the first kind relate 
the photon spectrum ("data") to the electron spectrum ("source function'). Such equa- 
tions are usually ill-posed in the sense of Hadamard - the effects of the experimental 
noise can b e strongly amplified by the intrinsic numerical instability of the model (for 
details, see iBertero et al.lll985l ). 

Standard approaches to obtain a soluti ons of inverse proble m in solar X-ray spec- 
troscopy are based on forward- fitting (e.g.. iHolman et al.ll20oj ') the photon flux spec- 
trum with parametric forms of the electron flux spectrum (see Section [4.1|) . However, 
in forward-fitting, the number of parameters utilized in the input form is generally 
small. This imposes severe, possibly artificial, constraints on the allowabl e form for the 
source function a nd is the main reason why inversion techniques (e.g., iBertero et al.l 
ll985l : IPianalll994l l. which find the best model-free non-parametric fit to the data sub- 
ject to physically sound constraints, are currently a very promising approach to data 
analysis in solar hard X-ray spectroscopy. 

A particularly promising technique is regularized inversion(e.g.. lB"ertero et al.lll985l l. 
The essence of the regularization technique is to seek a least-squares solution of the 
pertinent integral equation (e.g.. Equation |2.4| 1 within a subset of the solution space 
which accounts for some measure of a priori information on the source function. Con- 
sider the linear system 

g = Af, (4.2) 
which represents a discretized version of the (Volterra) integral equation (I2.4|l . where 

A, (^^^^, ) 5E,,^ = 1,...,N■ j-l,...,M; (4.3) 

f is the "source vector" nV{F{Ei), . . . , F(i5j\/)), g is the "data vector" {g{ei), . . . , g{ej^)) 
(with M > N), and the 5ei and SEj are appropriate weights. The values g{ei) corre- 
spond to a set of discrete photon counts in energy bands ej — >■ + Sei, while the F{Ej) 
are the corresponding values of the mean electron flux in energy bands Ej — !> Ej + SEj . 
Owing to the strong smoothing properties of the integral operators, the matrix A is 
quasi-singular and standard inver sion routines ca nnot be effectively applied. However, 
Tikhonov regularization theory (jTikhonovl [l963 ) obtains sufficiently smooth source 



functions as the (unique) solution of the minimization problem 

||g - Aff + A||Cf|p = minimum, (4.4) 

where A is the regularization parameter, which tunes the trade-off between the fit- 
ting term ||g — Af||^ an d the penalty term ||Cf|p. If C = I, the method is termed 



zero-order regularization ( Plana et al. l2003l l. while if C is the matrix cor responding to 



numerical differentiation, the method is termed first-order regularization ( Kontar et al.l 
[2004 V The optimal choice of the parameter A can be accomplished by means of some 
optimization approach or by means of a semi-heuristic technique based on a statistical 
analysis of the cumulative residuals ( Plana et al.ll200^ . 



The main disadvantage of using Tikhonov regularization is that for noisy data, 
solutions with negative (unphysical) values mi ght result. A possible so lution to this is 
provided by the projected Landweber method ( Plana fc Berter3ll997l ^ 



f„+i =P+[f„+TA^(g-Af„)] , fo = 0, (4.5) 
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where P+ sets to zero all negative components at each iteration and r is a relaxation 
parameter. In this framework, the tuning between stability and fitting is realized by 
applying some optimal stopping rule to the iterative procedure. 

4-2.1 Validation of regularization techniques 

The effectiveness of different in version algorithms, i ncluding a standard forward-fitting 
technique, have been tested bv lBrown et al. 1 l|2006l ) using synthetic data. In this test 



six forms of the mean source electron spectrum F{E) in Figure 14.31 (bottom), each 
one characterized by specific features like bumps or cutoffs, were used to generate the 
corresponding hard X-ray spectra 7(e) in Figure [53] (top), and a realistic amount of 
Poisson noise added. In all cases, the photon spectra look smooth and quite similar 
in their shape, while the corresponding mean electron spectra exhibit very irregular 
behavior that is filtered out by the smoothing effect of the bremsstrahlung cross-section. 

These differences epitomize the mathematical concept of ill-posedness. 

T he comparisons used four different techniques: zero- order (|Pianalll994l ') and first- 



order I Kontar et al.|[20o3 ) Ti khonov regularizatio n, triangular matrix row elimination 



with variable energy binning (jjohns fc Linlll992bh and forward-fitting with a paramet 



ric form consisting of a do uble power law with low- and high-energy cutoffs plus an 
isothermal component (e.g.. lHolman et al.ll2003l ). These te sts were done "in t he blind" 



to recover F(E) for later comparison with the true forms (|Brown et al.ll200d '). All of 
these approaches were able to reconstruct the general magnitude and form of F{E), 
although forward-fitting inevitably fails to recover small features which are not coded 
within the parameterized form of the model function. 

4.2.2 Application to RHESSI data 

The Tikhonov regularization method has been applied to hard X-ray measurements 

recorded by RHESSHot SOL 2002-02-20T11:07 (C 7.5), SOL2002-03-17T19:31 (M4.0), 

SOL2 002-08-06T12:59 (C7.9) dMassone et al.ll2003l). SOL20 02-02-26T10:27 (C9.6) l|Kontar et al.l 



2OO5I ) and SOL2002-07-23T00:35 fX4.8) (|Piana et al.ll2003l ). In this last paper a "dip 



hump" feature in the recovered mean source electron spectrum nVF was noted near 
E — 55 keV (Figure r4.4[l: such a feature is (by construction ) absent in t he superimposed 



forward-fit spectrum ( Holman et al.|[2003l ) using the same l|Hauelll997h bremsstrahlung 
cross-section. 



The 3(7 error bars plotted in Figure 14741 clearlv show that this "dip- hump" feature 
is statistically significant. Its physical interpretation is still an open issue, but it may 
reflect the depletion of low-energy nonthermal electrons due t o the effect of Coulom b 
collisions for an injected distribution with a low-energy cutoff (iEmslie fc Smithlll984l 'l. 
On the other hand, it may simply reflect an inadequate correction for the pulse pileup 
effect in the RHESSI detectors, w ith a resultant aliasing of the photon spectrum used 
in the construction of Figure SI] ijSmith et al.ll2002l ; iKontar et aLll2003h . 



4.3 High-energy cutoffs in the electron distribution 

It should be noted that the matrix A in Equation (|4.2p need not be square, so that 
the energy range corresponding to the electron flux (source) vector f may extend over 
a larger range than the photon flux (data) vector g. Physically, this corresponds to the 
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Fig. 4.3 Synthetic mean source electron spectra (bottom) and the corresponding photon 
spectra (top) in a "blind" experiment for assess ing the effectivenes s of different reconstruction 
methods in RHESSI X-ray spectroscopy (after iBrown et al.ll2006f) . 
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Fig. 4.4 Regularized spectrum nVF versus electron energy E for the time interval shown 
in SOL2002-07-23T00;35 (X4.8). The spectrum has been extended at high energies using a 
power law of index 5 = 2.45 (dashed lined). The vertical size of the error boxes reflects the 
3a- limit caused by statistical noise in the observed /(e). The spectrum obtained by a forward- 
fltting procedure using the same bremsstrahlung cross-section is shown as a solid line (after 
iPiana et aHl2003n . 
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production (through free-free bremsstrahlung only) of photons of energy e by electrons 
of energy E > emax, where emax is the largest photon energy observed. It should be 
noted that, because the bremsstrahlung cross-section Q{e, E) is a non-diagonal matrix 
of E, the form of the photon spectrum at photon energy e does provide information 
on the electron spectrum for all energies E > t, including those with E > emax- 
Consequently, a distinct advantage of the regulariz ation methodologies (over, say, the 
more straightforward matrix inversion method of I Johns fc Liiil Il992bl ') is that some 
information can be obtained on the form of the electron spectrum F{E) at energies 
beyond the maximum photon energy observed. 

A particular example of this is the possible existence of a high-energy cutoff -Emax 
in the electron spectrum F{E). Although for free-free emission this would correspond 
to a maximum photon energy emax ~ Smax, emax may lie b eyond the ran g e of s tatis- 
tically useful photon (or count) data. However, as shown bv lKontar et al.l (|2004h . the 
regularized inversion method can, in principle, detect the presence of the high-energy 
cutoff at iJ = -Emax through its effect on the photon spectrum at observable energies e 
that are all significantly less than -Emax, and this technique was indeed used to discern 
a high-energy cutoff (or at least a very sharp downw ard spectral break) in the fiare of 
SOL2002-02-26Tf0:27 (C9.6) bv lKontar et al.l (|2005'). The ability of the regularization 
technique to detect such high-energy c utoffs in the e l ectron spectrum was dramatically 
highlighted in the validation study of iBrown et al.l (|2006h . In that study, a synthetic 
electron spectrum F{E) with a high-energy cutoff at 200 keV (case B in Figure 14. 3p 
was used to generate noisy photon "data" in the range wholly below 100 keV. Anal- 
ysis of these "data" using both zero-order and higher-order regularization techniques 
rather faithfully reproduced the high-energy cutoff in the electron spectrum, with an 
accuracy better than 30%. High-energ y cutoffs have had very Umited use t o date as a 
parameter in forward- fitting methods ( Sui et a"l]|2002l : Holman et al.l 20031) . although 



both forward-fit and direct matrix inversion methods (|Johns fc Linlll992br) should be 
able to say something about the electron spectrum above maximum photon "data" 
energy, in that if there is a high energy cutoff -Emax close to the maximum observed 
photon energy the spectral shape is much different. 



4.4 Spectral breaks in the electron distribution 



Figure 1431 shows the photon spectrum for the time interval 09:43:16 - 09: 44:24 UT (the 
time of approximate peak flux) for SOL2005-01-17T09:52 (X3.8) (after iKontar et aD 
l2007l ). T his event, which produced several strong gamma-ray lines, was previously 
studied bv lKontar fc Brownl 1 2006bh . who concluded that the pitch-angle distribution 



for electrons up to ~300 keV was close to isotropic. There is evidence for an upward 
break in the spectrum at energies > 300 keV. For completeness, we note that the 
energy r ange > 300 also c ontains some pseudo-continuum 7-ray contributions from 
ions (e.g. ISmith et al.ll2003l ). 

The mean source electron spectrum F{E) recovered from the photon data, using 
two different bremsstrahlung cross-sections - electron-ion only and (electron-ion -f 
electron-electron) is shown in Figure 14.51 (The results are presented in the form of 
two "confidence strips," bundles of solutions using different noisy realizations of the 
same data.) The F{E) recovered using the full cross-section 1)2. 11[) . including electron- 
electron bremsstrahlung, is, for E >300 keV, steeper (spectral in dex greater by ^0.4) 
than the F{E) recovered assuming purely electron-ion emission. iKontar et al.l (|2007f ) 
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Fig. 4.5 Left panel: Photon spectrum for the time interval 09:43:16 - 09:44:24 UT in SOL2005- 
01-17T09:52 (X3.8), with gamma-ray lines removed. Right panel: Recovered forms of the quan- 
tity n V F{E) (in units of 10^" electrons cm~^ s~^ keV~^ ; see equation 12.41 ) using a zero-order 
regularization technique and presented as a "confidence strip," i.e., a series of solutions, each 
based on a realization of the data consistent with the size of the uncertainties. The dashed 
lines assume electro n-ion emission only; the solid lines include the additional electron-electron 
emission term (after iKontar et aHl2007r) . 



point out that while the upward break at e ~300 keV in Figure [431 is real, the break 
at iJ ~ 400 keV in the F{E) recovered using the electron-ion bremsstrahlung cross- 
section alone (Figure 14. 5p is an artifact of the neglect of electron-electron emission at 
energies ^300 keV. The true form of F[E), obtained using the full cross-section, has a 
rather straightforwa rd power -law form over the energy range from 200 — 1000 keV. As 
pointed out bv Kontar et al. 

1 f2007i ). the inclusion of electron-electron bremsstrahlung 
may remove th e need to explain "break energies" in seve ral events (cf. lDennislll985l : 
lHeristchilll986l : 019951 : iTrottet et aDll998l : IVestrandlll988l ') . However electron-electron 
bremsstrahlung cannot ex plain spectral breaks at low energies <300 keV often observed 
in RHESSI spectra (e.g.. (conwav et al.llioO^ '). 



4.5 Low-energy cutoffs in the electron distribution 

Given that nonthe rmal electron spectra often have a form close to a steep power law 
E~\ with 8 > 2 (|Dennislll985h . an accurate value for the low-energy cutoff Ec pa- 



rameter is required to obtain values of th e total energy in nonthermal electrons (see, 
however, lEmsliell2003l :l Hannah et al-lliooih . Thus, the determination of Eg plays a key 
role in the interpretation of hard X-ray data (see also lHolman et aDboill '). 

In most cases, the value of Ec must be somehow disentangled from the combined 
(thermal -I- nonthermal) form of F{E), and the value of Ec is, in general, only weakly 
constrained by observations. Some flares require rather high values of Ec to explain 
observations. The presence of a low-ener gy cutoff introduces a flattening of the photo n 
spectrum at ene rgies below Ec (see, e.g.. 



y cutolt mtroduces a nattenmg ot the photo n 
Holman|[2003 : lHuand l2009': Ha n et alEoogl l . 



Nitta et al.l ||199Q 1 show that the spectrum of the impulsive component flattens toward 



low energies, suggesting a value of Ec as high as 50 keV. iFarnik et al.l (|l997l ). using 
Yohkoh data, observed a few flares with flat spectra below 30 keV. 

Studies before RHESSI mostly assumed an arbitrary, flxed, value of Ec- More re- 
cently, using 7?_ff£S'S'/ high-spectral-resolution data (without albedo correction), several 
studies have reported clear evidence for a low-energy cutoff or even a dip (or gap) (e.g.. 
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iKontar fc Brownll2006al 'l in the mean electron flux distribution F{E) between its ther- 
mal and nonthermal components, leading to a clear identification of a low-energy cutofl' 
for the nonthermal electron distribution (Figure [4. 6|l . 

Forward-fitting methods (Section 14. 1|) usually assume a s trict low-energy cut off, 
i.e., a sharp change in the F{E). However, as showed by, e.g.. iKucera et a~ or 
iKasparova et al] (2003), it is often not possible to determine a unique best spectral fit; 
fit parameters can be varied substantially without unacceptably large changes in the 
photon spectrum, with other free parameters compensating. In particular, because of 
the dominance of thermal emission at low energies, forward-fit approaches can reliably 
infer only upper limits to the value of Ec; values of Ec well below such upper limits 
provide equally good fits to an observed X-ray spectrum, with a somewhat different 
value of the thermal source temperature T. Moreover, as mentioned in (Section I4.1|l 
determining confidence intervals for model parameters is an extremely time-consuming 
task. 

On the other hand, regularize d inversion metho ds (Section 14. 2[> have proven their 
ability to detect dips in F{E) (see lPiana eraDl2003l ). These methods also provide esti- 
mates of uncertainties in the solution, t hrough so-cal le d con fidence strips in which 
99% of data-consistent solutions lie; see iBrown et al. I (|200d ). Using this approach, 
clear low-energy cutoffs or dips in the me an electron fiux distrib ution F{E) at an 
energy around 20-40 keV were reported by iKasparova et al.l (|2005h for SOL2002-08- 
20T08 :25 (M3.4). F orward-fit methods yielded ^Bc = 44 ± 6 keV, a somewhat higher 
value. iKontar fc Brow n (2006a) found a clear dip around 20 keV for the SOL2002-09- 
17T05:54 (C2.0) event. For further examples of low-energy cutoffs obtained through 
both forward-fitting and regularization, see Figure 14.61 and Figure 14.71 respectively. 

Some flares, such as SOL2002-04-25T06:02 (C2.5), SOL2002-09-17T05:54 (C2.0), 
or SOL2002-08-20T08:25 (M3.4), have quite hard (flat) spectra with a re latively weak 
thermal components. Such a flat form of the photon spectrum can require ( Kontar fc Brownl 
[2006a) a low -energy cutoff , or local minimum, in the corresponding mean electron 
distribution F{E). These local minima are particularly interesting, since if F{E) is 
sufficiently steep (steeper than E^; see the remarks after Equation (|2.9|) '). it could have 
a form inconsiste nt with the widely-used collision-dominated thick-target model for 
X-ray production (i Kontar fc Brownll2006al : iBrow n ct al."2009''). 

However, the effects of Compton backscattering (Section I3.1|l o n the hard X-ray 
spectr um are most pronounced for flares with such hard spectra. Kasparova et al] 
(|2007l ) showed that the flares listed above are all located close to solar disk center 
{fj, > 0.5, 9 < 60°; denoted as stars in Figure and they therefore attributed the 
flat spectra to the heliospheric-angle -dependent albed o(Sec tion | 3.1l). This res ult is 
consistent with earlier observations o f iNitta et al ] (|l99d ) and lFarnfk et alJ (|l997l ). who 
observed several flares with flat spectra which were located not far from disk center 
(/X > 0.6) or from near center {fi m 1.0), respectively . Interestingly, adding considera- 
tions of the albedo in such events ( Kontar et al]|2()()6l ) removes the spectral hardening, 
and hence t he need for a low-energ y cut off in this photon en ergy range - see Figure 
14.61 and 14.71 IKasparova et aP ( 2005h , and IKontar et al ] (|2006h . Recently, IKontar et al] 
(|2008 a) have analyzed a large number of solar flares with weak thermal components 
and flat photon spectra. It has been shown that if the isotropic albedo correction is 
applied, all low-energy cutoffs in the mean electron spectrum are removed, and hence 
the low-energy cutoffs in the mean electron spectrum of solar flares above ~12 keV 
cannot be viewed as real features. If low-energy cutoffs exist in the mean electron spec- 
tra, their energies should be less than ~12 keV. Thus, the apparent low-energy cutoff 
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Fig. 4.6 Mean electron flux spectra nVF{E) of SOL2002-09-17T05:54 {C2.0) recovered using 
forwa rd fitting. The solid /dashed lines show the spectrum without /with the albedo correction 
(after iKontar et ani2006f) . 




10 100 



Energy, keV 



Fig. 4.7 Mean source electron flux spectrum nVF{E) for SOL2002-08-20T08:25 {M3.4) for 
the time interval 08:25:20 - 08:25:40 UT recovered using regularized inversion. The solid line 
shows the spectrum without albedo correction. The confldence intervals represen t the range of 
solut ions found for different statistical realizations of the photon spectrum (after lKontar et al.l 
l2006f) . 



in the mean electron distribution i s most likely to be a feature connected with albedo; 
it is not a true physical property ( Hannah et ahlEoogl l. This result can substantially 
change the total electron energy requirements in a given flare. 

In summary, our view on the existence and/or value of low-energy cutofTs has been 
significantly broadened since high quality RHESSI observations have become avail- 
able. Yet, the determination and even the existence of a low-energy cutoff remains 
a very complex issue. First, the contribution of the photospheric albedo to the ob- 
served photon spectrum must be taken into account and removed before the observed 
photon spectrum is converted to the electron distribution. Such correctio n removes ap- 
parent energy cutoffs in the mean electron distribution above ~12 keV ( Kontar et al.l 
bOOSal ). Secondly, apparent evidence from hard X-ray emission for a low-energy cut- 
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off in the electron spectrum should be carefully combined with other data sensitive 
to the low-ene rgy cutoff in nonthermal el ectron distribution, e.g., microwav e spectra 
( Holman|[20o3 ) or plasma radio emission (jMel'nik et all 1 19991 : iKontaillioOlh . Finally, 
effects which may lead to hard X-ray spectral flattening (e.g., albedo a nd non-uniform 
ionization; see lKontar et al.ll200l [iooa: iLiu et alll2009l : ISu et al.ll2009l ) should also be 
carefully assessed before reaching conclusions on the value of the crucial parameter Ec, 
with its attend ant implications for the energy content in nonthermal electrons (e.g., 
iFalewicz et"al][2009i V 



4.6 Temperature distribution of thermal plasma 

The thermal free- free (bremsstrahlung) continuum emission (photons s~^ keV~^) at 
photon energy e from an element of plasma of density n (cm~'^), temperature T (K) 
and volume dV is, neglecting factors of order unity, 

dI{e) = a^^eM~^/kT), (4-6) 

where k is Boltzmann's constant and a is a constant. For an extended source, with 
nonuniform density and temperature, integration of 1)4. 6|) over the spatial extent of the 
source gives 

^(^) = -J^ JTJ2 ^(2^) ^M-^/kT) dT, (4.7) 



where ^(T) = dV/dT (cm~ ^ K~^) is the differen tial emission measure at t emper- 
ature T (see equation (10) of ICraig &: Brownl Il976h . As first pointed out bv iBrownl 
(|l974l ). equation 1)4. 7|) may be written as a Laplace transform with respect to the 
inverse temperature variable x — 1/kT: 

— e7(6)=/ e'^'' f{x)dx = C[f{x);e], (4.8) 



where f{x) = ^{1/kx) x~^^'^ or, equivalently, ^(T) = (fcr)"^/^ /(l/fcT). 

A solution of equation ()4.8|l for f{x) (and hence (.(T)) is formally possible for a large 
variety of (but, it should be noted, not all; sec Brown & Emslie 1988) forms of /(e). 
Even when a solution does exist, however, the sol ution of equation (|4 .8I) is not a trivial 
task. This integral equation is of Fredholm-type dBertero et al ] |l985f r and is highly ill- 



posed, with a large class of solutions ^{T) corresponding to a given /(e) when (even very 
small) uncertainties in /(e) are taken into account. Mathematically, this extreme ill- 
posedness arises from the very broad form of the Laplace kernel exp(— ea;). Physically, 
the problem exists because of the broad range of temperatures T that contribute to 
the emission at a given photon energy e. Unlike, for example, the bremsstrahlung 
photon-to-electron inversion problem (see Section 14. 2p , in which only electrons with 
energy E > e contribute to the emission at photon energy e (and so the corresponding 
integral equation is of Volterra type, and so less ill-posed), in the electron-energy-to- 
temperature problem considered here all temperatures T contribute to the emission at 
energy e (and conversely a source at a single temperature T produces emission at all 
photon energies). 
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An early study of the inversion of equation (I4.8[l . us ing data from a ba lloon-borne 
instrument I Lin &: Schwartz! Ejirl ). was carried out by IPiana et all (|l995l l. They de- 



duced not only forms of ^(T) but also (in a one- dimensional geometry) the corre- 
sponding conduct ive flux -Fc(T) and its derivative dFc/dx. However, later work by 
IPrato et alJ ( 20061 ) examined the relation between (,{T) and the mean source electron 



spectrum F{E), viz., 

= -^m^M-E/kT)dT. (4.9) 

This study utilized the Landweber method, which ensures positivity of ^(r) everywhere 
- see equation |4.5) in Section 14.21 - and which is effective in recovering narrow, 5- 
function-like forms of the differential emission measure ^{T). It also showed that the 
form of ^(T) deduced from inversion of this (ri gorous) equa t ion is much less well- 
determined than suggested by the earlier work of IPiana et al.l l|l995l '). which involved 
inversion of the (inexact) Equation (|4.6|l . 

The method was applied to three photon spectra emitted during the flares SOL2002- 
08-21T01:41 (M1.4), SOL2003-11-03T09:55 (X3.9), and SOL2002-07-23T00:35 (X4.8). 
The photon spectra were first inverted by applying zero-order Tikhonov regulariza- 
tionto recover the form of F{E). For the first two events, the recovered ^{T) was 
consistent with a roughly isothermal low-temperature plasma plus a very broad form 
of C(r) at high temperatures. However, for the SOL2002-07-23T00:35 (X4.8) event, 
the reconstruction method produced unacceptably large residuals at low temperatures, 
consistent with the fac t that this same spectrum fails to satisfy the derivative test 
( Brown fc Emslielll98i) necessary for compatibility with a purely thermal interpreta- 



tion of the event. Observation of optically-thin lines is oft en used as an altern ative 
approach to infer the DEM of solar plasma (as a review, see IPhillips et aHl2008l ). 



5 The electron angular distribution 

Flare studies typically derive the properties of accelerated particles in the target region 
from observations of the radiation spectrum, but such radiation spectra are strongly 
dependent on the angular distribution of the energetic particles. Consequently, knowl- 
edge of both the angular distribution and the energy distribution of energetic particles, 
as they interact in the solar material, is necessary for understanding the acceleration 
and transport of particles in flaring regions. 

Attempts to measure the angular distribution of the accelerated electrons rely on 
the fact that an anisotropic ensemble of bremsstrahlung-producing electrons will gener- 
ate a radiation field that is both polarized and anisotropic. Efforts to measure electron 
beaming have therefore concentrated on studies of the hard X-ray continuum emis- 
sion, by looking at either the directivity or the polarization of the emitted radiation. 
Theoretical studies have considered the evolution of the electron pitc h-angle distri- 
bution as the part i cles are transported along magnetic field lines (e.g.. [Brownlll97"3 : 
iLeach fc PetrosianI Il98ll : iMcTiernan fc Petrosiaiil Il990al lbh. It is generally expected 
that, even if the particles are strongly beamed when injected, the net effect of the 
particle transport th rough the atmosphere will be to broaden the angular distribution 
(jHolman et al.ll201lh . Ideally, the measurements of the electron angular distribution 
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should therefore be performed as a function of both time and space. No such mea- 
surements are available as yet. However, those measurements that are available (as 
discussed below) already indicate some evidence for electron beaming. 

The angular distribution of accelerated ions can be studied by measu ring the en- 
ergies and widths of broad 7-ray lines (for details, see ISmith et al.ll2003l ). Studies of 
7-ray line data from the Solar Maximum Mission (SMM) Gamma Ray Spectrometer 
(GRS) suggest that proton s and a-particles are likely be ing accelerated in a rather 
broad angular distribution ( Share fc Murphvlll997l : IShare et al., ,2002). It is curren tly 
still unclear whether electrons and ions are being accelerated i Zharkova et ahlboilh in 
a similar fashion or by the same process. 



5.1 Early results 

One technique for studying hard X-ray directivity is to look for center-to-limb variations 
on a statistical basis. Correlations between flare longitude and flare intensity or spec- 
trum reflect the anisotropy of the X-ray emission and hence an associated anisotropy 
of the energetic electrons. For example, if the radiation is preferably emitted in a 
direction parallel to the surface of the Sun, then a flare located near the limb will look 
brighter than the same flare near the disk center. Various statistical studies of X-ray 
flares at energies below 300 keV report ed no significant ccntcr-to-limb variation of the 
observed intensity (e.g.. iDatlowe et al.l [l974: Piz zichini et al, 1974^ A statistically sig- 
nif icant center-to-limb va riation in the shape of the spectra of these events was found 
bv lRov fc Datlowel l|l975h . suggesting that perhaps some directivity may be present. 



The large sample of flares detected at energies greater than 300 keV by SMM /GRS 
allowed, for the first time, a statistical search for directivity at higher energies. Analysis 
of these data collected during Solar Cycle 21 provided the first clear evidence for di- 
rected emission, with a tendency for the high energy even ts to be located near the limb 



ected emission, witn a tendency tor tne nign energy even ts to be located near tne nmb 
Vestrand et ai] 1 19871 : [Bail 1 19881 ). iBoeovalov et al.l (|l985h and iMcTiernan fc PetrosianI 



199lh also reported evidence for anisotropies at hard X-ray energies. Obse rvations from 



5MM/GRS during Solar Cycle 22 provided further support for directivity (jVestrand et al.l 
Il99lh at energies above 300 keV. However, several high-energy events were also ob- 
served near the disk center by a number of different experiments during Cycle 22 (e.g., 
on GRAN AT an d the Compton Gamma-Ray Observatory, GGROh (For a summary, 
see IVilmeillTgiS ): this perhaps suggests a more complex pattern. |Lj (|l995l ) used data 
from SMM to confirm the general results of Vestrand et al. 

(1983) and concluded that 

there was evidence for increasing directivity with increasing energy. Quantifying the 
magnitude of the directivity from these statistical measurements is difficult. For exam- 
ple, one needs to know the size-frequency distribution for flares as w ell as the form o f 
the electron distribution to derive the predicted limb fraction (e.g.. |Petrosian|[l985l l. 
Furthermore, the results only represent an average for the flare sample. Different flar- 
ing regions are not likely to have identical geometries; nor are individual flares likely 
to have time-independent electron distributions. 

Another, more dire ct, method for studying dir ectivity in individual flares is the 
stereoscopic technique (jCatalano fc van Allenlll973l '). This method compares simulta- 



neous observations made on two spacecraft that view the same flare from different 
directions. ,Kane et a l. (1998) combined observations from the Pioneer Venus Orhiter 
(PVO) and ISEE-3 satellites to produce stereoscopic flare observations of 39 flares in 
the energy range from 100 keV to 1 MeV. While the range of flux ratios measured by 
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Fig. 5.1 Polar diagram of the bremsstrahlung cross-section lIGluckstern et al.lll953l ') for an 
electron of energy E = 100 keV and photon energies e = 30 keV (solid line), e = 50 keV 
(dotted line), and e = 80 keV (dashed line); the radial coordinate is proportional to the size of 
the cross-section and the angle from the z-axis corresponds to the angle between the incoming 
electron direction and the line to the observer. Note that at energies i? 2> e, t he cross-section 
peaks at = 0, while for i? ~ e, the cross-section peaks at 9 ~ (30 — 40)° (after lMassone et al.l 
12004) . 



Kane et al.l (|l998l ) is consistent with the results of statistical studies i Vestrand et al.l 
1987l l. the deviations of the ratio from unity show no clear correlation with increasing 



difference in viewing angles. Later stud ies concluded that there was no clear evidence 
for directivity at hard X-ray energies ( Kane et al.lll99i : iLi et al]|l994 ). Stereoscopic 
observations tend to suffer from cross-calibration issues between different instruments. 

The high quality hard X-ray data from RHESSI ha.ve opened new new opportunities 
and diagnostic techniques for the study of electron anisotropy in solar flares. 



5.2 Anisotropy of X-ray bremsstrahlung emission 

Using Equation H1.2|) . it is straightforward to show that the emitted bremsstrahlung in 
the direction ft toward the observer can be written 

/(e,f^) = ^^ j^T{E,n')Q{e,E,e')AEcm', (5.1) 

where Q{e,E,9') is the cross-section differential in photon energy e and the angle 6' 
between the precoUision electron direction fl' and the emit ted photon direction i7 , 
summed over the polarization state of the emitted photon taluckster n et al]|l953l l. 
The bremsstrahlung cross-section is more angle-dependent for higher photon energies 
(Figure [5. ip . which is qualitatively consistent with the findings of the statistical anal- 
ysis. 
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In Equation (|5.1|l the mean source electron flux spectrum has been generalized to 
the form F{E,fl'), which takes into account the angular distribution of electrons. If 
we denote by {9, (j>) the polar coordinates of the electron velocity vector ft' relative to 
the mean direction of the electron velocity distribution (usually the same direction as 
the guiding magnetic field line), and assume azimuthal symmetry, i.e., no dependence 
on (j), so that F{E, fj') = f{E, 6), and the angle of the mean direction of the electron 
velocity distribution is 9q, Equation (|5.ip can be written: 

— Tr I^OO fir flTT 

/(e) = -^/ / / Q[e,E,d'{e,<i,-do)\f{E,e)s.med^dedE , (5.2) 

where the directivity angle 6' between electron and photon directions is given by 
cos 6' = cos 9 cos 6q + sin 9 sin 9^ cos 4>. If we now define: 

/.27r 

Ko{e,E,9)^ Q[e,E,9'{9,^;9o)]d^ , (5.3) 

^0 = 

we can write 

^(^) = r^/ / Ko{e,E,9)f(E,9)sm9d9dE . (5.4) 

Inversion of Equation H5.4|) requires construction of the bivariate function f{E, 9) from 
knowledge of the (noisy) univariate function /(e). The problem is significantly more 
tractable if we make the further assumption that /(-E, 9) is separable in E and 9. We 
can then choose a particular form for the angular dependence of f{E, 9) and reconstruct 
the part of f{E,9) that depends only on the electron energy (univariate problem in 
E) or, analogously, we can assume the i?-dependence for f{E, 9) and recover the 9- 
depen dence (univaria t e pro blem in 9). 

In lMassone et all (|2004l l. the ^-dependence of f{E,9) is a prescribed (normalized) 
form g{9)/ J^, g{9) dfi'. Specifically, it was assumed that at all energies the pre-coUision 
electron velocities are uniformly distributed over a solid angle within a cone of half- 
angle a centered on a direction corresponding to a photon emission direction of 9q. 
With such an assumption. Equation 1)5. 4|) becomes 

/(e) = ^ 1^ (e, E) F{E) dE, (5.5) 

where ^ 

T^i^^E) = -—^ -/ Ko{e,E,9)sm9d9. (5.6) 

27r(l -cosa) Jg^^ 

Equation 1)5. 5p is formally identical to Equation (|2.4p . Application of the usual Tikhonov 
regularization technique to equation (|5.5|l then yields the mean source energy spectrum 

Fm , 

iMassone et al.' (iooi) applied this analysis to a photon spectrum recorded by the 
RHESSI instrument from SOL2002-08-21T01:41 (M1.4) in the time interval 00:39:04- 
00:39:48 UT (left panel of Figure [5TT1) . Values for the observation angle 9q ranging from 
0° to 180°, and for the cone semi-angle a ranging from 10° to 180°, were considered. The 
results (right panel of Figure 15. 2p demonstrated clearly that the angular dependence 
cannot be neglected: use of the anisotropic cross-section yields electron spectra that 
are significantly different from the ones reconstructed by using the solid- angle- averaged 
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e (keV) E (keV) 

Fig. 5.2 Left: photon flux spectrum. Right: regularized electron flux spectra corresponding 
to 00 = 130° (angle between an observer and the electron beam direction) and (from top to 
bottom) a = 10°, 30°, 60°, 90°, 180°. The choice of do = 130° is justified by the fact that 
this is the value corres ponding to a ver tically downward electron beam at the location of the 
selected flare (after iMassone et al.ii2004(X 



cross-section (corresponding to a = 180°). As the electron distribution becomes more 
anisotropic (decrease in a) , the cross-section for emission in the direction of the observer 
decreases for some parts of the cone and increases for others. The overall effect is a 
reduction in the cross-section, so that more electrons are required to produce the given 
photon spectrum. This effect is more pronounced at high energies (see Figure so 
the spectrum for small values of a is flatter than for a — 180° (isotropic distribution). 

As the viewing angle 8q increases (source moves closer to the disk center), the 
cross-section for emission in direction of the observer in general decreases, especially 
at high values of E, and accordingly the required electron spectrum increases and also 
flattens. 



5.3 Statistical results on X-ray anisotropy 

As a result of the heliocentric angle dependence of albedo (Section[3J , t he shape of pho- 



ton spectra should v ary as a function of their position on the solar dis k (iBai &: Ramatv 



19781 : iKontar et al] [2006). A statistical analysis of RHESSI flares (|KasDarova, et al 



20071 ) demonstrates a clear center-to-limb variation of photon spectral indices in the 
15 - 20 keV energy range and a weaker dependency in the 20 - 50 keV range. The 
observed spectral variations were found to be consistent with the predictions of albedo- 
induced spectral index changes (Figure [5. 3|l . 

Because the number of albedo photons depends on the amount of downward di- 
rected primary emission, the characteristics of the albedo component can be used to 
gt an esti mate of the a nisotropy of the primary hard X-ray emission. In this manner, 



get an esti mate o t the a nisotropy ot the primary hard A-ray emission. In this manner, 
iKasparova et al.l 1 2007h obtained values for a (the ratio of downward-directed primary 



emission to observer-directed primary emission) in the 15-20 keV energy range; see 
Figure [5741 The values of directivity a range from 0.3 to 3; while they imply downward- 
collimated emissio n in some cases, overall th ey are also cons istent with isotropic X-ray 



collimated emissio n in some cases, overall th ey are also cons istent with isotropic A-ray 
emission (see also iKontar fc Brownl [ioOGbl ) . The model of iLeach fc PetrosianI ( 19831 , 
Figure 4) predicts a < 3 at 22 keV for disk-center events. 
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Fig. 5.3 Spectral index 70 (15-20 keV) and 71 (20-35 keV) versus cosine fi of heliocentric 
angle 6. Vertical error bars indicate average uncertainties on the values as determined from 
single power-law fits. Lines show the predicted dependency for single power-law primary spectra 
with spectral indices 2.0, 2.5 and for directivity a{fi) = 1, 2,4 (solid, dashed, and dotted lines 
respectively). Flares with a dip in the mean electron distribution are denoted by stars; see also 
Section 14.51 (after iKasparova et al.ii2007l ) . 
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Fig. 5.4 Directivity at different heliocentric angle fi determined from the distributions of 
spectral index 70 in the 15 - 20 keV energy range at the limb and a given range of /i. Cross- 
hatched areas show the range of directivity values that are consistent (at the 0.05 significance 
level) with the hypothesis that the observed distribution at a given range of fi and the mod- 
eled distribution are drawn from the same parent distribution; crossed areas correspond to 
the 0.01 significance level. The dashed line shows the isotropic case, i.e., a{ij,) = 1 (after 
IKasparova et al.ll2007l') . 



5.4 Albedo as a probe of electron angular distribution 

The albedo spectral "contaminant" in fact offers very valuable insight into the anisotropy 
of the flare fast electron distribution. It does so by providing a view of the hard X-ray 
flare from behind, like a dentist's mirror. Moreover, the solar albedo "mirror" is spec- 
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trally distorting, so that its contribution to the overall spectrum can be distinguished. 
The observed spectrum in the observer's direction should contain an albedo "bump" 
feature, the strength of which is an indicator of the degree of downward beaming of the 
electron distribution. By use of this solar "mirror" we can achieve a degree of knowl- 
edge about the directionality of the primary photon distribution, and so the accelerated 
electron distribution, from single-spacecraft photon spectrometry ( Kontar fc Brownl 
l2006bh . 

The required integration over solid angle in Equation (|5.1|l can be approximated 
using a bi-directional representation. Invoking axial symmetry about a beam assumed 
to be perpendicular to the solar surface, we introduce the cross-sections 

cos(6io - AO) - cos(eo + A6) Jgg_^g 

averaged over [6q — A9,9q + A9] and centered at angle 9q. We define Q^{e,E) = 
Q(€,E,eo = 61) and Q^{e,E) = Q{e,E,9o = 180° - 61), where 9 is the heliocentric 
angle of the source. The electron spectrum F{E, 9) is defined in a similar bi-directional 
approximation, through introduction of the quantities 



nV J 



{Y)dV (5.8) 



for electrons propagating either upward (u) or downward (d) toward the scattering 
photosphere: Fu(E) and Fj^[E) are the density-weighted volumetric mean flux spectra 
of electrons directed towards the observer (upward and downward, respectively), also 
averaged over A9. 

With these assumptions and deflnitions, the discretized hard X-ray spectrum I(ej); i = 
1...A'', acco unting both for the prima ry spectrum lo and the reflected spectrum Ir, can 



be written (|Kontar fc Brownl [2006br ) 



I = Io+Ir = (Q* +G(/i)Q" qB+G(m)Q^) ( ) , (5.9) 

where F^, (j(i5j);j = 1...M axe the electron data vectors for downward- and upward- 
directed electrons, respectively and Q^'^ are matrix representations of the kernels of 
integral equations H5.7|) . The Green's matrices G(m) depend on the h e liocentric angle of 
the sou rce ^ = cos 9 and have been calculated in I Kontar et al.l ( 20061 ). iKontar fc Brownl 



(|2006b') have solved the inverse problem (Equation 15. 9|) using flrst-order Tikhonov 
regularization method. 

Figure [531 shows electron spectral solutions (Fu{E), F^liF)) for the M3.4 flare 
SOL2002-08-20T08:25. The results are consistent with isotropy up to £ 100 keV, 
with some indication of upward anisotropy above 100 keV. These flndings show a 
near-isotropic distribution (low ele ctron directivity) of the mean electron spectr um of 
accelerated electrons in solar flares (|Kontar fc Brown"2006bl: iHolman et al.ll201ll ). This 



strongly contradicts the models with pure l y coUisional transport in so l ar flares (e.g., 
BrownHl972l: IZharkova fc Gordovskvvl2006l : [Karlickv fc Kasparovall2009l : iHolman et al] 
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Fig. 5.5 Lower panel: The recovered mean electron flux spectra defined by Equation 15.81 
(thick lines) for SOL2002-08-20T08:25 (M3.4) (accumulation time interval 08:25:20-08:25:40 
UT). Downward-directed F^{E) (solid line) and observer- directed Fu{E) (dashed line) are 
shown, with corresponding 2cr errors (thin lines). Upper panel: electron anisotropy defined as 
Pd{E)/ Fu{E) with confidence values within the shaded area (after iKontar fc Brown. 2006bl) . 



5.5 X-ray polarization and electron angular distribution 



5. 5. 1 Model predictions 



The difficulties of statistical and stereoscopic observations for measuring hard X-ray 
directivity evoke the need for a technique that can measure anisotropics for individual 
flares. Because the bremsstrahlung cross-section Q{e,E,d') is in general polarization- 
dependent, polarization is a diagnostic that can, in principle, meet these requirements. 
Models of nonthermal (e.g., thick-target) hard X-ray production predict a clear and 
significant p olarization signal vyith polarization levels >10 % for beam- like distributions 
of electrons (iBai fc Ramatvlll978l:lBrownlll97ilCharikov et al.lll996l: lEmslie fc Vlahod 



1980l : iLaneer fc Petrosianll 19771 : iLeach fc Petrosianlll983l : IZharkova et al.lll995l ) and an 



orientation parallel to the plane containing the guiding magnetic field direction and the 
direction to the observer. For vertical guiding fields, this ori e ntation direction projects 



onto a radial line on the solar disk ( Bai fc Ramatvl 1 19781 : iLeach fc PetrosianI Il983l : 



IZharkova et al.lll995h : for other orientations of th e guiding field, other polarization 
vector orientations are possible jEmsIie et al.l l2008ll. For strong electron beaming, the 



polarization degree can r each values up to 60 % at energies above 50 keV, or even 
higher for lower energies (|Hauelll972l ). 
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All models predict a strong dependence between the observed value of polarization 
and the viewing angle. The highest polarization degrees are expected for large angles of 
view, when the line of sight is perpendicular to the magnetic field line. Thus, most the- 
ories predict higher polarization for flares located near the solar limb. The direction of 
the polarization depends on the energy: low energy ha rd X-rays are negativ ely polar- 



ized whereas, above 350 keV, the sign turns to positive i Bai fc Ramatvl[l978l ). Because 
collisions of beam electrons with ambient particles tend to isotropize the distribution, 
the highest polarizations (up to 85%) are expected at ene rgies around 100 keV ge n- 



erated in the coronal portions of the top of the flare loop ( Leach fc PetrosiarJIlQSj l. 
Photons observed from the footpoints (in the region of the dense chromosphere) would 
be polarized only to the level of around 20%. It should also be noted that even thermal 
models of the hard X-ray source predi ct a finite polarization of order a few percent, 
due to thermal gradients in the source l|Emslie fc Brownlll98ol '). The thermal compo- 



nent, with its rather low polarization, tends to dominate the emission from all flares 
at energies below about 25 keV. 

All these theoretical predictions, while clearly testable, could be criticized on the 
grounds that the modeling assumptions they contain may be oversimplistic. For ex- 
ample, each model to date assumes a single, simple magnetic field structure about 
which the emitting electrons spiral. It could be argued that any real flare, particularly 
one sufficiently large to produce a signal of sufficient strength to enable a polarization 
measurement, will in all probability contain a mix of structures that would average 
out any polarization signal present. However, hard X-ray imaging observations in the 
impulsive phase generally show a fairly simp le geometry, consisting of two footpoin t 
sources and perhaps a loop-top source (e.g., iMasuda et aPllQQSl : ISakao et al.lll992l '). 
These observations suggest that simple magnetic structures are responsible for the 
energetic emissions and give support to the possibility that a statistically significant 
polarization signal could be produced in a large event. 

As noted in Section [3l and in other places throughout this article, a substantial 
fraction of the observed hard X-ray flux is backscattered from the solar photosphere, 
as the so-called photo sphericalbedo. The precise magnitude of this backscattered frac- 
tion depends, in part, on the polarization of the primary flux. Further, the reflected 
component will influence the degree of polarization of the total observed fiux, since 
backscattering will t end to introduce polarization fractions of a few percent at energies 
below 100 keV (e.g.. lBai fc Ramatvl[l978l : iLanger fc Petrosianll 19771 ). Direct imaging of 
this albedo patch would place a constraint on the contribution of such backscattered 
photons to the primary sig nal. iHudson et all (|2003l ) have also suggested that Compton 
scattering in the corona may lead to measurable polarization effects. Clearly, a simul- 
taneous hard X-ray imaging capability (such as that provided by RHESSI) represents 
a major advantage for interpretation of a hard X-ray polarization measurement. 

5.5.2 Htstory of observations 

The history of observations of hard X-ray polarization from solar ffares is a fascinating 
subject in its own. The first measurements of X-ray polarization from solar ffares (at 
energies of ~15 keV) were made by Soviet experimenters using polarimeters aboard 
the Intercosmos satellites. Their polarimeters were made of a hexagonal Be-scatterer 
surrounded by six counters located in front of the Be-prism faces. Later versions of the 
instrument were mounted on a turnable drum to reduce systematical errors. In their 
initial study, [rindo et aL 1 lll970h reported an average polarization of P = 40% (±20%) 
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for three 1969 X-ray flares: SOL1969-10-20T10:50 (C9.0), SOL1969-10-23T05:15, and 
SOL1969-10-30T09:30 (M7.3). This study was followed by an analysis of two flares 
in 197 0: SOL1970-10-24T05:41 (M6.4) and SOL1970-11-05T03:21 (X2.3) (|Tindo et alJ 
Il972al lbh that showed polarizations of approximately 20% during the impulsive phase. 
These reports were met with considerable skepticism, on the grounds that they did not 
adequately allow fo r proper detector cross-calibration and had limited photon statistics 
( Brown et al]|l974l ). Subsequent observations with an instrument on the 050-7 satel- 
lite seemed to confirm the existence and magnitudes o f the polarizations ( ^10%), but 
these data were compromised by in-flight gain shifts jNakada et al.lfl97i ). In a later 
study using a polarimeter on Intercosmos 11, iTindo et'liLr(fl976ll measured polariza- 
tions of only a few percent at ~15 keV for two flares in July 1974. This small but 
finite polarization is consistent with the predictions for pur ely thermal emission that 
contains an admixture of polarized backscattered radiation ( Bai fc Ramatvl[l978l ). 

A decade later, a new polarimeter, designed to measure in the energy range from 5 
keV to 20 k eV (with about 1.5 keV energy resolution), was fiown on the Space Shuttle 
Columbia ( Tramiel et al.|[l984h . This design was based on metallic lithium scattering 
elements surrounded by gas-proportional counters. Contamination of the Li scattering 
elements invalidated the pre-fiight instrument calibration. An in-fiight calibration was 
performed using two flares near the center of the Sun, under the assumption of null 
polarization. The upper limits to the polarization derived from this calibration for 6 
flares (C and M-classes) were in the range from 2.5% and 12.7% (99% confldence). 

Recent measurements at energies below 100 k eV have been perfo rmed with the 
SPR-N instrument on the CORONAS-F satellite (jZhitnik et aklbood l. The SPR-N 
instrument included a solar X-ray radiation monitor and a polarimeter capable of 
detecting signals in the energy ranges of 20-40, 40-60 and 60-100 keV. The polarization 
detector consisted of a hexahedral Be-scatterer and three pairs of CsI(Na) scintillation 
detectors located on the faces of the Be prism. With a total effective area ranging 
between ~0.3 cm^ at 20 keV to ~1.5 cm^ at 100 keV, it detected hard X-rays from 
more than 90 flares between 2001 and 2005. From a sample of 25 solar flares, one 
could determine the upper limits of the polarization degree from 8 to 40% (3cr). Only 
for the single case of SOL2003-10-29T20:49 (XIO.O), located near Sun center, was a 
significant polarization level measured: its value increased from 50% (20-40 keV) to 
greater than 70% (60 — 100 keV). Although the same flare was observed by RHESSI, 
an independent polarization analysis was not possible because RHESSI was located 
at high magnetic latitude and was experiencing a high level of charged particle events 
(fSuarez- Garcia et al..i2006 ). 



5.5.3 RHESSI polarization measurements 

RHESSI enables polarization measurements in a wide range of energies fro m 20 keV to 



1 Me V. Its detection system of 9 cylindrical coaxial germanium detectors ([Smith et al.l 
|2002| ). coupled with a satellite rotation every 4 seconds, strongly reduces systematic ef- 
fects in the polarization measurements. However, the RHESSI design is not optimized 
for studies of polarization and its small effective area for this kind of study, together 
with the high background contribution, have led to measurements with limited statis- 
tical signiflcance. 

For the energy range between 20 keV and 100 keV, polarization can be measured 
using the photons that are scattered into the rear segments of the Ge detectors by 
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Fig. 5.6 Solar flare polarization data from lBoggs et al] ll2006l ) for photon energies between 
200 keV and 1 MeV. The best fit curve is shown as a solid line. The case for 100% polarization 
is shown as a dotted line. The measured polarization for SOL2002-07-23T00:35 (X4.8) is 21% 
± 9%. The measured polarization for SOL2003-10-28T11;10 (X17.2) is -11% ± 5%. 



a (passive) Be scattering block that is pla ced between four of th e germanium detec- 
tors within RHESSFs spectrometer array jMcConnell et al.ll2002l ). At energies above 
~100 keV, polarization measurements are performed using the so-called "coincidence 
mode": photons that Compton-scatter from one Ge detector to another are identified 
by a suitable coincidence timing window. In this way, one detector plays the role of an 
active scattering element, significantly reducing the background level. 

iMcConnell et all (|2003l ) reported on a measurement of polarization covering the 
energy range from 20 keV up to 100 keV performed with the Be scattering block. 
The initial results from an analysis of SOL2002-07-23T00:35 (X4.8) showed evidence 
for polarization at a level of 15±4%. The measured polarization angle (as measured 
counterclockwise from solar west) was 'I' — 79° (±5°), implying that the polarization 
vector is inclined ~64° with respect to the radial direction at the flare site. While the 
magnitude of the observed polarization is broadly consistent with the prediction of solar 
flare models which invoke the precipitation of a nonthermal electron beam into a dense 
chromospher ic target, the orientation of the polarization vector is somewhat surprising. 
lEmslie et al.l ( 2008 ) suggest that the orientation of the polarization vector in this case 
could be explained by a tilt of the flaring loop with respect to the local v ertical. Such 
a tilt is also consistent with gamma-ray line observation s for this flare ('Smith e t al.l 
l2003l l. Unfortunately, recent analysis of these data (Mc Connell et a l. 2007) suggests 
that at least a part of the reported signal may be a result of systematic effects in the 
data analysis, underscoring the difficulty of making reliable polarization measurements. 
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Using the RHESSI coincidence mode, two results have been recently published 
for energies above 100 keV. The first one (see Figure 15. 6|) describes the polarization 
measurements of two X-class solar fla res, one located clo se to the solar limb and the 
other close to the center of the disk ( Boggs et al.11200^ ). The energy range selected 
was between 200 keV and 1 MeV. The polarization degrees found were 21% ± 10% 
and —11% ± 5% in la respectively. Assuming that the measured polarization was 
significant in each case, they measured a radial polarization direction for the flare close 
to the solar center, and an azimuthal direction for the fiare near the limb. The levels of 
polariza tion, as well as their di rections, are consistent with beamed electron distribution 
models ('Bai fc Ramatv"l978^) and contrary to the results using the albedo method 
(fkontar & Brown 2006b; K asparova et al.l 20071 ). However, the level of inconsistency 



with more isotropic distributions is rather weak. Given that a report of GRB (gamma- 
ray burst) polarization by this same team has been met with considerable skepticism in 
the literature, one might also question these fiare results. It should be noted, however, 
that the fiare analysis differed significantly from the GRB analysis and has not been 
contested in the literature. 



The second study (see Figure \5?7} applies the coincidence met hod to the impulsive 
phase of seven solar flares in the 100 keV to 350 keV energy range I Suarez-Garcia et al.l 
I2OO6 ). The fiare sample consisted of six X-class (from XI. 4 to X8.4) flares and one 
M7.0-class flare located either on the limb or in the outer part of the solar disc. Val- 
ues for the polarization degree were 2% and 54%. The lowest degree of polarization 
obtained was found to be fully compatible with a 0% polarization measurement. The 
highest degree of polarization, which was the most statistically significant of the seven 
results obtained, was found to be 2.6 sigma away from 0. The angles of polarization 
'I' were distributed between 35° and 85° independen t of the fiare location, contrary 
to both the results mentioned in iBoggs et al.l 1 20061 ). and to the expectation based 
on simple geometrical modeling. Additional attempts to correlate various parameters 
(e.g., polarization level, polarization angle, heliocentric angle, footpoint orientation, 
flare intensity) were also i nconclus ive. The results were com pared with the theoretical 
predictions from lBai fc Ramatv (19781 ) and lLeach fc Petrosian (1983) and with the 0% 
polarization hypothesis. Th e analysis allowed only for rejection (90% of confldence) 
of one of the models from iLeach fc PetrosianI (|l983l ). In this model, predicting very 
high polarization values up to 85%, the magnetic field strength is constant along the 
loop and the electrons spiral at pitch angles close to 90°. Due to the statistical un- 
certainties, for the rest of the models the values were very close to unity, making 
it im possible to d i stingu ish b etween them. The polarizat ion amplitudes measured by 
both lSoggs et all ([2006) and ISuarez-Garcia et al.l (|2006l ) are combined in Figure \E 
with diamonds and filled circles, respectively. 



The results to date do not yet provide unambiguous evidence for solar flare polar- 
ization at hard X-ray energies. It is clear, however, that such data would constitute a 
unique probe into the electron acceleration process in solar flares, providing new con- 
straints on theoretical models and allowing for more detailed studies of the acceleration 
processes and geometries. 
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Fig. 5.7 Solar flare polarization data from [Suarez-Garcia et al\ ll2006fl for photon energies 
between 100 and 350 keV. The best fit curve is shown as a solid line. The best case for 
polarization is from SOL2005-01-19T08:22 (X1.3), where a polarization value of 54%±21% 
was measured. 
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Fig. 5.8 Results on the absolute value of the degree of hard X-ray polarization above 100 keV 
achieved with RHESSI (Icr errors). The diamonds and the filled circles correspond to two inde- 
pendent measurements (see text for references). From low to high flare class (left to right) the 
flares represented are: SOL2005-08-25T04:40 (M6.4), SOL2005-01-19T08:22 (XI. 3), SOL2004- 
11-10T02:13 (X2.5), SOL2005-01-17T09:52 (X3.8), SOL2002-07-23T00:35 (X4.8) (two mea- 
surements), SOL2005-01-20T07;01 (X7.1), SOL2003-11-02T17:25 (X8.3), and SOL2003-10- 
28T11:10 (X17.2). None of these results are significant at a level of greater than 3cr, suggesting 
that higher quality data are still needed. 



6 The electron spatial distribution 

Hard X-ray imaging spectroscopy is a powerful tool with which to explore the underly- 
ing physics of particle acceleration and transport in solar flares, and has been a central 
component of the RHESSI concept since its beginnings. 



6.1 Early results 

In its most basic form, imaging spectroscopy involves simply constructing and com- 
paring two-dimensional (count) maps of the source for different energy bands. The 
earliest imaging of solar hard X-ray sources was carried out using the Hard X-Ray 
Imaging Spectrometer (HXIS) (|van Be ek e t al.i,198Q ) on (SMM). This instrument pro- 
vided imaging information between 3.5 and 30 keV though the use of an array of 
subcoUimator-defining grids; it had an angular resolution of 8 arcsec over a 1 60 arc - 
sec field of view, and 32 arcsec over a wider, 444 arcsec field. iHovng et al.1 (|l981al ) 



showed that the hard X-ray emission in SOL1980-04-07T01:07 (M4) was located "in 
two patches;" the patch with the harder spectrum "coincided with the brightest Ha 
emission." On the other hand, in SOL1980-04-10T09:23 (M4), the hard X-ray emis- 
sion was "concentrated in a looplike structure, with the softer spectrum at the top of 
the loop and the harder spectrum in the legs, thus ind i cating preference for the thick- 



target model of hard X-ray production." iHovng et al.l 1 1981bl ) studied the hard X-ray 
emission from a large two-ribbon flare SOL1980-05-21T20:50 (XI) and concluded that 
the higher-energy (16-30 keV) emission originated in "separate locations of ~8 arcsec 
width, coinciding in position with Ha flare kernels," while the softer (3.5 — 8 keV) 
emission originated from a "broader region in between." 
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TreatinR these o bservations as evidence for a liot coronal region ( Krucker et alJ 
l2008l : lTomczai3l2009l l produced by the primary energy release process, and el ectron pre- 
cipita tion from this primary acceleration region into the dense chromosphere. lMachado et alJ 
ijigsd ) studied the energetics of SOL1980-04-10T09:23 (M4) and conclud ed that "only a 
fractio n of the total flare energy" was present in the accelerated electrons. IPuiiveman et ahl 
(|l982l ') continued the study of this event, in addition to SOL1980-05-21T20:50 (XI) and 
SOL1980-11-15T15:53 (XI). By comparing the energy in the accelerated electrons with 
an estimate of the thermal conductive flux out of the hot coronal acceleration region, 
they concluded that "a large fraction of the dissipated flare power has to go into elec- 
tron acceleration." 

The Solar X-ray Telescope (SXT) instrument I Takakura et al.lll983l ) on the Hino- 
tori satellite used the rotating modulation colli mator (RMC) te chnique with two pairs 
of rotating grids to make hard X-ray im a ges JOdaet al.lll976h: th e FWH M angular 
resolution was ~30 arcsec. [Tsuneta et al.l (| 19831 ) and Tsuneta et al.l ( 19841 ) report ob- 



servations of a near-limb flare, SOL1981-05-13T06:10 (XI. 5), which ex hibited a diffuse 
hard X-ray source situated some 40,000 km above the photosphere. iTakakura et al.l 
( 19861 ) reported that a hard X-ray (20 — 40 keV range) image for the impulsive com- 
ponent of SOL1982-02-22T04:44 (M2.7) wa s an extended source elongated along the 
solar limb with a source height of 7,000 km. iTsuneta et al.l (Il983l ) also report imaging 
hard X-ray observations of SOL1981-07-20T14:41 (M5.5), for which "each hard X-ray 
source in the initial phase coincides with each Ha bright region," further evidence for 
a thick-target interpretation of the hard X-ray emission at such energies. 

In contrast to the RMC technique use d by Hinotori and R HESSI, the Ifard X-ray 
Telescope (HXT) instrument on Yohkoh ( Kosugi et al.|[l99ll ) used information from 
64 pairs of occultation grids; each pair of grids yields a single F ourier component 
of the source. The FWHM resolution achievable was ~8 arcsec. [Matsushita et all 
(|l992l ). using observations of about a hundred flares observed with the Yb/ifco/i/HXT, 
reported that, on average, the hard X-ray source height "decreased with increasing 
X-ray energy," c onsistent w ith the deeper pen etration of higher-energy electrons (e.g., 
iBrown fc McClvmont, 197ll : lEmslielll978l . ll98ll ). With nine RMC components, an angu- 
lar resolution d own to ^2 arcse c, and high spectral resolution, the information available 
from RHESSI (|Lin et al.ll20o3 ) is far superior for imaging spectroscopy studies. 



6.2 Imaging spectroscopy with RHESSI 

One o f the earlie st results using RHESS I imeigmg spectroscopy data was bv lAschwanden et al.l 
(I2OO2I ): see also iKrucker fc LinI (|2002l ). who analyzed hard X-ray sour ce height as a 
function of energy, improving the previous stati stical analvsis (lM atsusliit a~et al.lll99^ ) . 
With the superior data available from RHESSI, lAschwanden et al.l (|200 j ) were able to 
determine accurate source locations as a function of photon energy for the single event 
SQL200 2-02-20Tll:07 (C7.5), ra ther than a statistical average for an ensemble of events 
used bv [Matsushita et al ] l|l992l ). This was done by determining the centroid location 



of a circular Gaussian forward-fit to each of the two clearly resolved footpoints in this 
event. Then, using thick-target modeling of the expected centroid location vs. pho- 
ton energy in a planc-parallel-stratified atmosphere model I Brown et al.ll20o3 : iMroze^ 



I2OO6 : Kontar ot al. 2008b), they were able to deduce a density vs. height structure 
for the atmosphere that was consistent with empirical models of flaring atmosphere 
density. 
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As mentioned at the outset of this section, imaging spectroscopy is, in principle. 



at different energies (e.g., Kruckcr & Lin 2003: Saint-Hilairc & Bena 2002: Sui et al. 
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2007; 
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Shao & Huand I2OO9I: Petrosian & Chen I2OI0I). These maos 


are produced by apply- 



ing image processing algorithms (e.g., back -projection, CLEAN, Maximum Entropy or 



Pixon; for details see iHurford et aLll2002l ) to the temporally modulated fluxes from 
each of RHESSFs detectors, in which spatial information on the source is encoded. 
Then, "interesting" regions in the field of view are selected and the intensity in such re- 
gions is determined in the map corresponding to each energy range. There results from 
this a set of intensity-versus-energy profiles, i.e., a count spectrum, for each feature. 
Using knowledge of the instrument response permits photon spectra for each feature to 
be determined. Finally, the corresponding spatially-resolved electron spectra are con- 
struc ted through forwar d fitting or by applying regularized spectral inversion methods 
(e.g.. lBrown et al]|2006l ) to the spatially resolved photon spectra. Analysis of the vari- 
ation of the electron spectrum throughout the target is a powerful dia gnostic of the 
phys ical processes affecting the bremsstrahlung-producing electrons (see lEmslie et all 
I2OOII ). 



Emslie et al.l 



(|2003l l performed such a "stacked-image" analysis for different energy 
bands via an imaging spectroscopy analysis of SOL2002-07-23T00:35 (X4.8). They iden- 
tified four "interesting" features in the images (Figure 16. ip and constructed a count 
spectrum for each (Figure 16. 2p . They identified one of these features as a "coronal" 
source, and showed that its spectrum was consistent with a thermal source at a temper- 
ature of order 45 x 10^ K. The spectra of the other three sources were more power-law- 
like Jn_Jorm^_riuSj_£lus their location near the lower altitudes in the flaring structure, 
1120031 ) to identify these features as chromospheric footpoints. 

"South" 



ledlEmslie et al.l 



Interestingly, the spectra of two of these footpoints (labeled "North" and 
in Figure [S7T|, while significantly different and v arying with time, m aintained a rela- 
tively constant difference throughout the event. lEmslie et al.l 1 2003h interpreted the 
similar time variation to a magnetic coupling of the two footpoints, and the system- 
atic difference in spectral indices as due to differential spectral hardening associated 
with an asymmetric location of the electron acceleration region. From the magnitude 
(~0.3) of the spectral-index difference, they obtained an estimate of the differential 
column density between the acceler a tion r egion and each of the footpoints . More re- 
cently, however, ISaint-Hilaire et aP l|2008l ). in a statistical study of some 50 events, 
have pointed out that such a large intervening column density should in general (but, 
somewhat ironically, not necessarily in SOL2002-07-23T00:35) produce more emission 
in the legs of the loop than was observed by RHESSI; they conclude that electron pre- 
cipitation in asymmetric magnetic field geometries is a more reasonable explanation 
for the observed footpoint spectral differences. 

While the above imaging spectroscopy technique appears relatively straightforward 
to implement, there are a num ber of dif ficulties t hat warrant some commentary. First, 
for large events, pulse pileup ( Datlowg 1976b,a: iSmith et al.|[2002l ') is an issue. Pairs 
of low energy photons arriving nearly simultaneously are detected as a single energy 
count. Since the modulation of such counts corresponds to the low energy source, an 
image made at the higher energy returns a "ghost" image of the low energy source. 
Second, the dynamic range of RHESSI is such that features containing a few percent 
of the total flux cannot be reliably imaged, so that spectral information in the relevant 



Deducing electron properties 



47 




9 "^--^ j 60-120 keV 

0040 



36-41 keV 



northern 
middle tm, 



41-46 keV 



75-85 keV 



96 109 keV 



□ 



109 123 keV 1S3 139 keV 139-156 keV 156-177 keV 




Fig. 6.1 Top panel: time profile for event. Lower panels: Images of SOL2002-07-23T00:35 
{X4.8) in different count energy channels, for time interval 3 of the event (cf., flux proflle in 
top panel). The four main regions in the event are labeled. Since this e vent was located n ear 
the East limb, structures to the left are higher in the atmosphere (after lEmslie et al.ll200li) . 



energy ranges is not reliably recovered. For example, in the lEmslie et al.l (l2003h analysis 
of SOL2002-07-23T00:35, the spectrum of the "coronal" source was determined only 
at energies < 40 keV, while the spectrum of the footpoints was determined only at 
energies > 30 keV. Only in the relatively narrow energy range from 30 to 40 keV could 
the spectra of all the features in the source be obtained. 



6.3 Visibilities and imaging spectroscopy 

It is important to realize that, because of the RMC technique used by RHESSI, spatial 
information is encoded in the RHESSI data in a distinctive way, namel y in rapid time 



variations of the detected counts in each of the RHESSI subcoUimators. Schmahl et al.l 
have developed a technique in which the observed temporal modulations pro- 
duced are interpreted in terms of a set of vistbiliUes (calibrated measurements of specific 
spatial Fourier components of the source distribution). As with image reconstruction 
in radio interferometry, the set of visibilities thus determined can then be used to infer 
the spatial properties of the X-ray source. It is important to note (1) that visibili- 
ties are a "first-order" product of the RHESSI data, and (2) that the data can be 
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Fig. 6.2 Spectra for each of the features identified in Figure 16.11 For each time interval, 
the "coronal" source is well-fit to a thermal spectrum with the temperatures shown, and the 
"footpoint" sources are well-fit by power laws with the spectral indices shown. The spectral 
indices of the various footpoints vary with time; however, the spectral inde x difference betwee n 
the "North" and "South" footpoints is relatively constant with time (after lEmslie et al.ll200dl . 



used to determine not only the values of the visibilities, but also their quantitative 
uncertainties. 

IXu et al have analyzed a set of ten events, each of which exhibits a rather 

simple, single-extended-source, geometry. They performed their imaging spectroscopy 
analysis by assuming a parametric form of the source structure (a seven-parameter 
curved elliptical Gaussian) and then forward fit not to the actual images themselves, 
but rather to the corresponding source, visibilities Because both the visibilities and 
their uncertainties were determined quantitatively, they were able to deduce not only 
the best-fit valu es of the source param eters, but also the uncertainties on the value of 
each parameter. iKontar et al.l ( 2008bh have performed such an analysis with circular 
Gaussian fits (four parameters) of footpoints for the SOL2004-01-06T06:29 (M5.8) fiare 
and have measured the sizes and heights of the hard X-ray sources in a few energy 
ranges between 18 and 250 keV. The height variations of footpoint emission with 
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energy and the source size with height have been found with an unpreced ented vertical 
resolution of ~150 km a t chromospheric heights of 400 — 1500 km (cf., iBrown et alj 
l2002l : [Aschwanden et al.| [2002'). T he interpretation of these findings in terms of electron 
transport models is given bv lHolman et al] 1 201ll ). 

One of the major goals of hard X-ray imaging spectroscopy is the determination 
not just of the source structure as a function of photon (or count) energy, and hence 
the variation of the hard X-ray spectrum throughout the source, but rather of the vari- 
ation of the corresponding electron spectrum throughout the source. It is the spectral 
variation of this electron spectrum that provides the key insight into the physics of 
electron acceleration and propaga tion in solar flares. 

Recognizing this, IPiana et al.l ([2007) introduced a new approach to imaging spec- 
troscopy which is optimized to the distinctive way in which spatial information is 
encoded in the RHESSI data. Although it is possible to use the count visibilities deter- 
mined rather straightforwardly from the raw data to construct count images and then 
procee d to do imag i ng spe ctroscopy analysis in the "stacked image" manner described 
above, IPiana et al. (|2007l ) point out that such "traditional" imaging algorithms are 
completely ineffective in smoothing in the energy direction, with the result that recov- 
ered images corresponding to adjacent energy bins can exhibit substantial differences. 
Further, RHESSFs Fourier-component approach to imaging detects "patterns" of emis- 
sion, rather than information in a "pixel-by-pixel" format, so that analysis of a partic- 
ular sub-region of a so urce is affect e d by the signal (and noise) contained in all other 
features in the source. iPiana et al. therefore argue that imaging spectroscopy 

analysis using RHESSI data is best accomplished through conversion of the temporal 
modulations in terms of count visibilities and subsequent analysis of these count visi- 
bilities to obtain information on the electron spectrum in the spatial-frequency (rather 
than spatial) domain. This leads to a set of electron visibilities, which contain all the 
information on the variation of the electron spectrum throughout the source (albeit as 
a function of spatial frequency, rather than position). If desired, such electron visibil- 
ities can then be used to construct electron flux maps using the same algorithms used 
to convert count visibilities into count maps. 

Figure 16.31 illustrates the essence of this procedure. Visibilities in count space are 
used to determine the corresponding set of electron visibilities using the T ikhonov 
regularization technique (Tikhonov 1963i) previously used (e.s.. IpTana et ali r2003) to 
ascertain spatially-integrated electron spectra F{E) from observations of spatially- 
integrated count (or photon) spectra 7(e) (see Section [4. 2p . Applied to visibilities, the 
Tikhonov regularization method permits the determination of information on ele ctro n 
fluxes at energies above the maximum count energy observed ( Kontar et al.l l2004l ) . 
More importantly, it forces smoothness in the inferred electron visibility spectra at each 
point in the spatial frequency domain and so enhances real features that persist over a 
relatively wide energy band, while suppressing noise-related features that show up only 
over a narrow range of energies. The combination of visibility data and the Tikhonov 
regularization methodology allows the derivation of the most robust information on 
the spatial structure of the electron flux spectrum image, the key quantity of physical 
int erest. 

IPiana et al] (|2007l ') have illustrated the power of the electron visibility method 
by applying it to data obtained near the peak of SOL2002-02-20T11:07 (C7.5). Us- 
ing visibilities from RHESSI RMCs 3 through 9, corresponding to spatial resolutions 
from ~7" to ~183", they construct the amplitude and phase of the count visibilities 
V{u,v;q) (counts cm~^ keV~^) as a function of Fourier components {u,v) and 
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Fig. 6.3 Construction of electron flux images from measured visibilities. Panel (a) shows the 
"traditional" approach to image reconstruction, in which the visibility information in count 
space from each of RHESSFs nine detectors is used to construct an image based on hard X-ray 
counts. Note that there is no imposition of energy smoothing in this process: images in adjacent 
count energy channels (q^ ± Aqi) can exhibit significant differences due to count statistics. 
These differences would be further magnified if the count-ba sed images were in verted to yield 
electron images. Panel (b) shows the innovative approach of iPiana et al ] 1I2OO7I I. in which the 
energy spectrum fo r each count visib ility (left panel) is subjected to a regularized inversion 
procedure (see, e.g.. lPiana et ani2003ri to yield the corresponding (smooth) electron spectrum 
for each visibility. Once all the electron visibilities have been determined (right panel), they 
can be used to yield electron images using the same image reconstruction algorithms as used 
to produce count-based images (Panel (c)). These electron flux images are, by construction, 
necessarily smooth across electron energy ranges (Ej ± AEj), and so are more amenable to 
further analysis. 
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Fig. 6.4 Electron flux spectral images corr esponding to the regularized electron flux spectral 
visibilities, obtain ed through application of iBong et alj 's (2006) MEM-NJIT algorithm (after 
iPiana et aLll2007D . 



count energy q. They then apply the Tikhonov inversion algorithm to construct the 
corresponding set of electron visibilities W{u,v; E). By construction, these electron 
visibilities vary smoothly with electron energy E. 



Figure [63] shows electron flux images for the illustrative SOL2002-02-20T11:07 
(C7.5) event. They show evidence for two footpoints, connected by a "strand" of coronal 
flux. Because of the inherent smoothness demanded by the regularized algorithm used 
to construct electron visibilities from co unt visibilitie s , the electron flux images vary 
smoothly with energy. As pointed out bv lPiana et al. I (|2007l ). this contrasts markedly 
with the behavior in the count images, for which the image in each energy range is inde- 
pendent, so that statistical fluctuations, including those in sidelobes from neighboring 
features, result in a set of images that do not vary smoothly with energy. Inversion of 
the count spectra obtained from such images leads to ampliflcation of such noise and 
electron spectra that exhibit large (and most probably) unphysical features. 



Three different spatial subregions in the source are highlighted in Figure 1631 Two 
of these regions correspond to the footpoint sources visible at higher energies and 
the other one to similarly-sized regions located approximately midway between the 
two footpoints. The lower panel of Figure [S3] shows the areally-averageqj electron-flux 
spectra (electrons cm~^ keV~^ arcsec"^), for each of these three subregions. These 
spectra are sufficiently smooth that significant conclusions regarding the variation of 
the electron spectrum throughout the source can be made. 



^ To get the total count spectrum for each region [counts cm ^ s ^ keV simply multiply 
the areally-averaged spectrum by the area of that region, viz., 14.4 X 14.4 = 207.36 arcsec-^ 
(Footpoint 1), 22.8x9 = 205.2 arcsec^ (Middle), and 14.4x14.4 = 207.36 arcsec^ (Footpoint 2), 
respectively. 
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Fig. 6.5 Top panels: Electron images in the energy ranges 22 — 26 keV and 42 — 46 keV, 
respectively. Three sub-regions of interest are labeled on each image. Two of these correspond 
to bright footpoint-like sources and one to a region midway between the footpoints. Bottom 
panel: Areally- averaged electron fl ux spectra (electro ns cm"'^ s~^ keV~^ arcsec"^) for each of 
the three sub-regions shown (after [Piana et aLll2007fl . 



7 Summary 

7.1 Hard X-Ray emission processes 

With the launch of RHESSI, not only have new types of hard X-ray data analysis 
become possible, but also the in-depth study of X-ray producing processes has been 
triggered (Section [2]). The roles of electron-ion and electron-electron bremsstrahlung, 
free-bound electron-ion emission, and Compton backscatter of primary photons, have 
been highly scrutinized in view of the unprecedented quality of the RHESSI data. 
Examples include: 

• As pointed out in Section [2.31 pure electron- ion bremsstrahlung spectra have a 
spectral index 7 ~ 5 -f 1, for pure electron-electron bremsstrahlung the resulting X- 
ray spectrum has a significantly shallower photon spectrum, with 7 ~ 5. Hence, the 
importance of the electron-electron bremsstrahlung contribution increases with photon 
energy and the enhanced emission per electron leads to a flattening of the photon 
spectrum /(e) produced by a given F{E), or, equivalently, a steepening of the F{E) 
form required to produce a given /(e). 
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• The recent work on the importance of long-neglected free-bound emission by 
energetic electrons (Section l2.2[l . with its emphasis on emission due to electron capture 
onto various ionization states of various elements, necessitates further work on the 
ionization structure of the solar atmosphere and its time dependence during flares. 

• The solar atmosphere above an X-ray emitting region can be safely treated as an 
optically-thin medium, whereas the lower levels of the atmosphere are optically thick 
for X-rays. X-rays at energies below ~11 keV are mostly photoelectrically absorbed 
in the photosphere, while Compton scattering dominates at the energies above. As a 
result. X-rays emitted downwards can be Compton-backscattered toward the observer; 
as discussed in Sections 13. II and 13. 21 they play a major role in the spectral and imaging 
characteristics of the emission. 



7.2 Electron source spectrum 

The spatially-integrated photon spectrum is a key source of information about the 
mean electron flux spectrum in solar flares. High-energy-resolution data from RHESSI 
have permitted, for the first time, not only the extraction of basic parameters through 
forward-fitting the observed spectra (Section 14. but have also permitted the first 
reliable model-independent inversion of observed spectra (Section 14. 2p . Although the 
general form of a large number of nonthermal fiare spectra can be adequately approx- 
imated by an isothermal Maxwellian for the low-energy component (Section 14. 6p . plus 
a broken power law for the high-energy component, RHESSI d&ta, have clearly demon- 
strated a wealth of features beyond such simple parametric models: the presence of 
high-energy and low-energy cutoffs (Sections l4.3l and l4.5|l . spectral breaks (Section l4.4p . 
and the presence of an albedo component leading to "dips" in the mean electron spec- 
trum (Section I4.2|l . The inclusion of an isotropic albedo correction removes the need 
for low-energy cutoffs and, if low-energy cutoffs exist in the mean electron spectrum, 
they should be below ~12 keV (Sections 13.11 and I4.5|l . 



7.3 Anisotropy 

The electron-ion bremsstrahlung emission cross-section is generally anisotropic, with a 
dipole-like diagram at low energies (Section (5]). The anisotropy or angular distribution 
of X-ray emitting electrons can be measured in a number of different ways, and var- 
ious possibilities have been employed with RHESSI data. Parameter-free regularized 
electron flux spectra, reconstructed by assuming a parameterized form of the electron 
angular distribution, highlight the need to consider anisotropy in determining the true 
shape of the electron flux spectrum (Section 15. 4|l . A statistical study of center-to-limb 
variations in the 15-50 keV energy range has shown that an anisotropy factor a, or 
the ratio of downward and upward directed fluxes, for hard X-ray emission that lies 
outside the range [0.2, 5] can be rejected with 99% confldence (Section 15. 3|) . 

The albedo portion of the observed spectrum can be rather effectively used to infer 
the mean electron flux in two directions simultaneously (Section l5.4|l . The reconstruc- 
tions of the mean electron spectra in the downward and upward directions suggest 
that X-ray emitting electrons need not be signiflcantly anisotropic in a broad range of 
energies from tens of keV to about 200 keV. This imposes signiflcant challenges to solar 
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flare m odels based on coUisional transport; see lFletcher et aP (|201lh and lHolman et aTl 
(|201lh for details. 

RHESSI has been used to attempt polarization measurements in two different en- 
ergy bands using two different scattering processes. The results are, however, of lim- 
ited statistical significance. If real, they would suggest that X-ray polarization above 
50 keV is significant (with values up to 50% or more) and that high quality polariza- 
tion measurements would provide significant constraints on particle acceleration models 
(Section I5.5|l . The RHESSI design is not optimized for studies of polarization, and 
attempts to measure it have been hampered by the small effective area and the high 
background contribution. Future progress will require a polarimetric instrument able 
to measure polarization with errors on the level of 1 to 2%. Such measurements would 
allow detailed studies of the electron beaming, quantifying its magnitude, and placing 
significant constraints on the acceleration geometry. 

7.4 Spatial variation of electron fiux 

Imaging spectroscopy results (Section[6ll from RHESSI have (not surprisingly) revealed 
that the hard X-ray emission from solar flares is far from homogeneous: the centroid 
of high energy emission moves downward with higher energy, presumably due to the 
increased penetration of higher energy electrons; coronal sources have radically different 
spectra than chromospheric footpoints; and even systematic differences exist between 
different footpoint sources in the same event. 

Since RHESSI provides spatial information on detected hard X-ray emission through 
a (Fourier-transform-based) rotating modulation collimator technique, it follows that 
the highest-fidelity spatial information is contained in the finite number of Fourier com- 
ponents ("visibilities") sampled. Both forward-fit and inversion techniques have been 
applied to RHESSI imaging spectroscopy data, with very positive results, including 
the first empirical estimates of the density and volume of the electron acceleration re- 
gion, and the recovery of electron flux maps that, by construction, vary smoothly with 
energy and so provide valuable information on the variation of the emitting electron 
spectrum throughout the fiare source. 
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